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ABSTRACT 

In the history of astronomy, major advances in our understanding of the Universe frequently arise 
from dramatic improvements in our abiUty to accurately measure astronomical quantities. Aided by 
rapid progress in information technology, current sky surveys are changing the way we view and study 
the Universe. Next-generation surveys will maintain this revolutionary progress. We focus here on 
the most ambitious survey currently planned in the visible band, the Large Synoptic Survey Telescope 
(LSST). LSST will have unique survey capability in the faint time domain. The LSST design is driven 
by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar 
System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a large, wide- 
field ground-based system designed to obtain multiple images covering the sky that is visible from 
Cerro Pachon in Northern Chile. The current baseline design, with an 8.4m (6.5m effective) primary 
mirror, a 9.6 deg^ field of view, and a 3.2 Gigapixel camera, will allow about 10,000 square degrees 
of sky to be covered using pairs of 15-second exposures in two photometric bands every three nights 
on average, with typical Scr depth for point sources of r ^ 24.5. The system is designed to yield high 
image quality as well as superb astrometric and photometric accuracy. The survey area will include 
30,000 deg^ with 5 < -f34.5°, and will be imaged multiple times in six bands, ugrizy, covering the 
wavelength range 320-1050 nm. The project is scheduled to have first light in 2014 and the beginning 
of survey operations in 2015. About 90% of the observing time will be devoted to a deep-wide-fast 
survey mode which will observe a 20,000 deg^ region about 1000 times (summed over all six bands) 
during the anticipated 10 years of operations, and yield a coadded map to r ~ 27.5. These data 
will result in databases including 10 billion galaxies and a similar number of stars, and will serve the 
majority of science programs. The remaining 10% of the observing time will be allocated to special 
programs such as a Very Deep and Fast time domain survey. We illustrate how the LSST science 
drivers led to these choices of system parameters, and describe the expected data products and their 
characteristics. The LSST data products will be available to the public and scientists around the 
world - everyone will be able to view and study a high-definition color movie of the deep Universe. 

Subject headings: astronomical data bases: atlases, catalogs, surveys — Solar System — stars — the 
Galaxy — galaxies — cosmology 
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1. INTRODUCTION 

Major advances in our understanding of the Universe 
have historically arisen from dramatic improvements in 
our ability to "see". We have developed progressively 
larger telescopes over the past century, allowing us to 
peer further into space, and further back in time. With 
the development of advanced instrumentation - imag- 
ing, spectroscopic, and polarimetric ~ we have been able 
to parse radiation detected from distant sources over 
the full electromagnetic spectrum in increasingly subtle 
ways. These data have provided the detailed information 
needed to construct physical models of planets, stars, 
galaxies, quasars, and larger structures. 

Until recently, most astronomical investigations have 
focused on small samples of cosmic sources or individual 
objects. This is because our largest telescope facilities 
typically had rather small fields of view, and those with 
large fields of view could not detect very faint sources. 
With all of our existing telescope facilities, we have still 
surveyed only a minute volume of the observable Uni- 
verse. 

Over the past two decades, however, advances in tech- 
nology have made it possible to move beyond the tra- 
ditional observational paradigm and to undertake large- 
scale sky surveys. As vividly demonstrated by surveys 
such as the Sloan Digital Sky Survey (SDSS; York et al. 
2000), the Two Micron AU Sky Survey (2MASS; Skrut- 
skie et al. 2006), and the Galaxy Evolution Explorer 
(GALEX; Martin et al. 2006), to name but a few, sensi- 
tive and accurate multi-color surveys over a large fraction 
of the sky enable an extremely broad range of new scien- 
tific investigations. These projects, based on a synergy of 
advances in telescope construction, detectors, and above 
all, information technology, have dramatically impacted 
nearly all fields of astronomy - and many areas of fun- 
damental physics. In addition, the recent world-wide at- 
tention received by Google Sky^'^ (Scranton et al. 2007) 
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demonstrates that the impact of sky surveys extends far 
beyond fundamental science progress and reaches all of 
society. 

Motivated by the evident scientific progress enabled by 
large sky surveys, three nationally-endorsed reports by 
the U.S. National Academy of Sciences'^'* concluded that 
a dedicated ground-based wide-field imaging telescope 
with an effective aperture of 6-8 meters is a high pri- 
ority for planetary science, astronomy, and physics over 
the next decade. The Large Synoptic Survey Telescope 
(LSST) described here is such a system. The LSST will 
be a large, wide-field ground-based telescope designed to 
obtain multi-band images over a substantial fraction of 
the sky every few nights. The survey will yield contigu- 
ous overlapping imaging of over half the sky in six optical 
bands, with each sky location visited about 1000 times 
over 10 years. 

The purpose of this paper is to provide an overall sum- 
mary of the main LSST science drivers and how they led 
to the current system design parameters (§ [2]), to de- 
scribe anticipated data products (§ [3]), and to provide 
a few examples of the science programs that LSST will 
enable (§11]). The community involvement is discussed in 
§ [51 and broad educational and societal impacts in § [51 
The concluding remarks are presented in § [7l This publi- 
cation will be maintained at the arxiv.org site^^, and will 
also be available from the LSST website (www.lsst.org). 
The latest version should be consulted for the most up- 
to-date information about the LSST system. 

2. FROM SCIENCE DRIVERS TO REFERENCE DESIGN 

The most important characteristic that determines the 
speed at which a system can survey a given sky area to a 
given depth (faint flux limit) is its etendue (or grasp), the 
product of its primary mirror area and the field-of-view 
area (assuming that observing conditions such as seeing, 
sky brightness, etc., are fixed). The effective etendue 
for LSST will be greater than 300 m^ deg^, which is 
more than an order of magnitude larger than that of 
any existing facility. For example, the SDSS, with its 
2.5-m telescope (Gunn et al. 2006) and a camera with 
30 imaging GGDs (Gunn et al. 1998), has an effective 
etendue of only 5.9 m^ deg^. 

The range of scientific investigations which will be en- 
abled by such a dramatic improvement in survey capabil- 
ity is extremely broad. Guided by the community-wide 
input assembled in the report of the Science Working 
Group of the LSST^^, the LSST is designed to achieve 
goals set by four main science themes: 

1. Probing Dark Energy and Dark Matter 

2. Taking an Inventory of the Solar System 

3. Exploring the Transient Optical Sky 

4. Mapping the Milky Way 

3* Astronomy and Astrophysics in the New Millennium, NAS 
2001; Connecting Quarks with the Cosmos: Eleven Science Ques- 
tions for the New Century, NAS 2003; New Frontiers in the Solar 
System: An Integrated Exploration Strategy, NAS 2003. 

35 http://lanl. arxiv.org/abs/0805.2366 

3S Available as |http://www.lsst.org/Science/docs/DRM2.pdf| 
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Each of these four themes itself encompasses a vari- 
ety of analyses, with varying sensitivity to instrumental 
and system parameters. These themes fully exercise the 
technical capabilities of the system, such as photomet- 
ric and astrometric accuracy and image quality. About 
90% of the observing time will be devoted to a deep- 
wide-fast (main) survey mode. The working paradigm 
is that all scientific investigations will utilize a common 
database constructed from an optimized observing pro- 
gram (the main survey mode), such as that discussed in 
Section 3. Here we briefly describe the science goals and 
the most challenging requirements for the telescope and 
instrument that are derived from those goals, which will 
inform the overall system design decisions discussed be- 
low. For a more detailed discussion, we refer the reader 
to the LSST Science Requirements Document'^^, as well 
as to the numerous LSST poster presentations at the re- 
cent meetings of the AAS'^*. 

2.1. The Main Science Drivers 

The main science drivers are used to optimize various 
system parameters. Ultimately, in this high-dimensional 
parameter space, there is a manifold defined by the to- 
tal project cost. The science drivers must both justify 
this cost, as well as provide guidance on how to optimize 
various parameters while staying on the cost manifold. 

Here we summarize the dozen most important inter- 
locking constraints on data properties placed by the four 
main science themes: 

1. The depth of a single visit (an observation con- 
sisting of two back-to-back exposures of the same 
region of sky) 

2. Image quality 

3. Photometric Accuracy 

4. Astrometric Accuracy 

5. Optimal exposure time 

6. The filter complement 

7. The distribution of revisit times (i.e., the cadence 
of observations), including the survey lifetime 

8. The total number of visits to a given area of sky 

9. The coadded survey depth 

10. The distribution of visits on the sky, and the total 
sky coverage 

11. The distribution of visits per filter 

12. Data processing and data access (e.g., time de- 
lay for reporting transient sources and the software 
contribution to measurement errors) 

We present a detailed discussion of how these science- 
driven data properties are transformed to system param- 
eters below. 

^'^ Av ailable at | htt pY ^www.lsst.org^Scieiice/docs.shtml| 

38 See |http://www.lsst.org/Meetings/AAS/2008/AAS211.shtml| 



2.1.1. Probing Dark Energy and Dark Matter 

Current models of cosmology require the existence of 
both dark matter and dark energy to match observa- 
tional constraints (Spergel et al. 2007; Riess et al. 2007; 
Tegmark et al. 2006; and references therein). Dark en- 
ergy affects the cosmic history of both the Hubble ex- 
pansion and mass clustering. If combined, different types 
of probes of the expansion history and structure history 
can lead to tight constraints (a percent level precision) 
on the dark energy equation of state and other cosmo- 
logical parameters. These tight constraints arise because 
each technique depends on the cosmological parameters 
or errors in different ways. The most powerful probes in- 
clude weak gravitational lens cosmic shear (WL), baryon 
acoustic oscillations (BAO), supernovae (SN), and clus- 
ter counting - all as functions of redshift. Using the cos- 
mic microwave background as normalization, the com- 
bination of these probes can yield the needed precision 
to distinguish among models of dark energy (Zhan 2006, 
and references therein). In addition, time-resolved strong 
galaxy and cluster lensing probes the physics of dark 
matter. This is because the positions and profiles of mul- 
tiple images of a source galaxy depend sensitively on the 
total mass distribution, including the dark matter, in the 
lensing object. 

Three major programs from this science theme pro- 
vide unique and independent constraints on the system 
design: 

• Weak lensing of galaxies, 

• Baryon Acoustic Oscillations 

• Type la Supernovae. 

Weak lensing (WL) techniques can be used to map the 
distribution of mass as a function of redshift and thereby 
trace the history of both the expansion of the Universe 
and the growth of structure (e.g., Hu & Tegmark 1999; 
Wittman et al. 2000; for a review see Bartelmann & 
Schneider 2001). Weak lensing of galaxies as a function 
of redshift probes both the evolution of structure over 
cosmic time and determines dimensionless ratios of dis- 
tances versus cosmic time, giving multiple independent 
constraints on the nature of dark energy. These investi- 
gations require deep wide-area multi-color imaging with 
stringent requirements on shear systematics in at least 
two bands, and excellent photometry in all bands (a re- 
quirement shared with BAO.) The strongest constraints 
on the LSST image quality come from this science pro- 
gram. In order to control systematic errors in shear mea- 
surement, it is mandatory to obtain the desired depth 
with many short exposures (which enables reconstruction 
of galaxy shapes with minimum systematic error). De- 
tailed simulations of weak lensing techniques show that, 
in order to obtain a sample of ^3 billion lensing galaxies, 
the coadded map must cover ~20,000 deg^, and reach a 
depth of r ~ 27.5 {5a for point sources), with several 
hundred exposures per field and sufficient signal-to-noise 
in at least five other bands to obtain accurate photomet- 
ric redshifts. This depth optimizes the number of mea- 
sured galaxies in ground-based seeing, and allows their 
detection in significant numbers to beyond a redshift of 
two. It is anticipated that optimal science analysis of 
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weak lensing will place strong constraints on data pro- 
cessing software, such as simultaneous analysis of all the 
available images rather than analyzing a deep coadded 
image (Tyson et al. 2008). 

Type la supernovae (SN) provided the first evidence 
that the expansion of the Universe is accelerating (Riess 
et al. 1998; Perlmutter et al. 1999). To fully exploit 
the supernovae science potential, light curves sampled in 
multiple bands every few days over the course of a few 
months are required. This is essential to search for sys- 
tematic differences in supernova populations (e.g., due 
to progenitor channels) which may masquerade as cos- 
mological effects, as well as to determine photometric 
redshifts from the supernovae themselves. Unlike other 
cosmological probes, even a single object can provide 
useful constraints and, therefore, a large number of SN 
across the sky can enable a high angular resolution search 
for any dependence of dark energy properties on direc- 
tion, which would be an indicator of new physics. Given 
the expected SN flux distribution, the single visit depth 
should be at least r ~ 24. Good image quality is re- 
quired to separate SN photometrically from their host 
galaxies. Observations in at least five photometric bands 
are necessary to ensure that, for any given supernova, 
light curves in several bands will be obtained (due to the 
spread in redshift). The importance of K-corrections to 
supernova cosmology implies that the calibration of the 
relative offsets in photometric zero points between filters 
and the knowledge of the system response functions, es- 
pecially near the edges of bandpasses, must be accurate 
to about 1% (Wood-Vasey et al. 2007). Similar photo- 
metric accuracy is required for photometric redshifts of 
galaxies. Deeper data (r > 26) for small areas of the 
sky can extend the discovery of SN to a mean redshift of 
0.7 (from 0.5 for the main survey), with some objects 
beyond z ^1 (Garnavich et al. 2005; Pinto et al. 2005). 
The added statistical leverage on the "pre-acceleration" 
era would improve constraints on the properties of dark 
energy as a function of redshift. 

2.1.2. Taking an Inventory of the Solar System 

The small-body populations in the Solar System, 
such as asteroids, trans-Neptunian objects (TNOs) and 
comets, are remnants of its early assembly. The his- 
tory of accretion, coUisional grinding, and perturbation 
by existing and vanished giant planets is preserved in the 
orbital elements and size distributions of those objects. 
Collisions in the main asteroid belt between Mars and 
Jupiter still occur, and occasionally objects are ejected 
on orbits that may take them on a collision course with 
Earth. 

As a result, the Earth orbits within a swarm of as- 
teroids; some number of these objects will ultimately 
strike the Earth's surface. In December 2005, the U.S. 
Congress directed'^^ NASA to implement a near-Earth 
object (NEC) survey that would catalog 90% of NEOs 
larger than 140 meters by 2020. About 20% of NEOs, 
the potentially hazardous asteroids or PHAs, are in or- 
bits that pass sufficiently close to Earth's orbit, to within 
0.05 AU, that perturbations with time scales of a cen- 
tury can lead to intersections and the possibility of col- 
lision. In order to fulfill the Congressional mandate us- 

For details see |http: //neo.jpl.nasa.gov/neo/report2007.html| 



ing a ground-based facility, a 10-meter class telescope 
equipped with a multi-gigapixel camera, and a sophis- 
ticated and robust data processing system are required 
(Ivezic et al. 2007a). The search for NEOs also places 
strong constraints on the cadence of observations, requir- 
ing closely spaced pairs of observations (two or preferably 
three times per lunation) in order to link observations 
unambiguously and derive orbits. Individual exposures 
should be shorter than about 30 seconds each to minimize 
the effects of trailing for the majority of moving objects. 
The images must be well sampled to enable accurate as- 
trometry, with absolute accuracy of at least 0.1 arcsec. 
The images should reach a depth of at least ~24.5 (5cr for 
point sources) in the r band in order to probe the 0.1 
km size range at main-belt distances, and to fulfill the 
Congressional NEO mandate. The photometry should 
be better than 1-2% to enable color-based taxonomy. 



2.1.3. Exploring the Transient Optical Sky 

Recent surveys have shown the power of measuring 
variability for studying gravitational lensing, searching 
for supernovae, determining the physical properties of 
gamma-ray burst sources, probing the structure of active 
galactic nuclei, and many other subjects at the forefront 
of astrophysics (Tyson 2006, and references therein). A 
wide-area dense temporal coverage to deep limiting mag- 
nitudes would enable the discovery and analysis of rare 
and exotic objects such as neutron star and black hole 
binaries, novae and stellar flares, gamma-ray bursts and 
X-ray flashes, active galactic nuclei (AGNs), stellar dis- 
ruptions by black holes, and possibly new classes of tran- 
sients, such as binary mergers of black holes (Shields & 
Bonning 2008). Such a survey likely would detect numer- 
ous microlensing events in the Local Group and perhaps 
beyond, and open the possibility of discovering planets 
(e.g., Beaulieu et al. 2006), as well as obtaining spectra 
of lensed stars in distant galaxies. 

Time-domain science requires large area coverage to 
enhance the probability of detecting rare events; good 
time sampling, because light curves are necessary to dis- 
tinguish certain types of variables and in some cases to 
infer their properties (e.g., determining the intrinsic lu- 
minosity of Type la supernovae depends on measure- 
ments of their rate of decline); accurate color informa- 
tion to assist with the classification of variable objects; 
good image quality to enable differencing of images, espe- 
cially in crowded fields; and rapid data reduction, clas- 
sification and reporting to the community in order to 
flag interesting objects for spectroscopic and other in- 
vestigations with separate facilities. Time scales ranging 
from 1 min, to constrain the properties of fast faint tran- 
sients (such as optical flashes associated with gamma-ray 
bursts; Bloom et al. 2008) and transients recently dis- 
covered by the Deep Lens Survey (Becker et al. 2004), 
to 10 years to study long-period variables and quasars 
(Kaspi et al. 2007) should be probed over a signiflcant 
fraction of the sky. It should be possible to measure col- 
ors of fast transients, and to reach faint magnitude limits 
in individual visits (at least the Deep Lens Survey limit 
of r ~ 24.5). 
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2.1.4. Mapping the Milky Way 

A major objective of modern astrophysics is to under- 
stand when and how galaxies formed and evolved. Ar- 
guably the biggest challenge to extragalactic cosmology 
today concerns the formation of structure on sub-galactic 
scales, where baryon physics becomes important, and 
where the nature of dark matter may manifest itself in 
observable ways. The Milky Way and its environment 
provide a unique data set for understanding the detailed 
processes that shape galaxy formation and for testing the 
small-scale predictions of our standard cosmology. How- 
ever, we still lack robust answers to two basic questions 
about our Galaxy: 

• What is the detailed structure and accretion his- 
tory of the Milky Way? 

• What are the fundamental properties of all the 
stars within 300 pc of the Sun? 

Key requirements for mapping the Galaxy are large 
area coverage, excellent image quality to maximize 
the photometric and astrometric accuracy, especially in 
crowded fields; photometric precision of at least 1% to 
separate main sequence and giant stars (e.g., Helmi et 
al, 2003); and astrometric precision of about 10 mas 
per observation to enable parallax and proper motion 
measurements. In order to probe the halo out to its 
presumed edge at ^100 kpc (Ivezic et al. 2003) using 
numerous main-sequence stars, the total coadded depth 
must reach r > 27, with a similar depth in the g band. To 
study the metallicity distribution of stars in the Sgr tidal 
stream (e.g., see Majewski et al. 2003) and other halo 
substructures at distances beyond the presumed bound- 
ary between inner and outer halo (~30 kpc, Carollo et 
al. 2007), the coadded depth in the u band must reach 
~ 24.5. To detect RR Lyrae stars beyond the Galaxy's 
tidal radius at ~300 kpc, the single-visit depth must be 
r ^ 24.5. In order to measure the tangential velocity of 
stars at a distance of 10 kpc, where the halo dominates 
over the disk, to within 10 km/s needed to be competi- 
tive with large-scale radial velocity surveys, the required 
proper motion accuracy is at least 0.2 mas/yr. The same 
accuracy follows from the requirement to obtain the same 
proper motion accuracy as Gaia (Ferryman et al. 2001) 
at its faint limit (r 20) . In order to produce a complete 
sample of solar neighborhood stars out to a distance of 
300 pc (the thin disk scale height), with 3(t or better ge- 
ometric distances, trigonometric parallax measurements 
accurate to 1 mas are required. To achieve the required 
proper motion and parallax accuracy with an assumed 
astrometric accuracy of 10 mas per observation per co- 
ordinate, approximately 1,000 observations are required. 
This requirement on the number of observations is in 
good agreement with the independent constraint implied 
by the difference between the total depth and the single- 
visit depth. 

2.1.5. A Summary and Synthesis of Science-driven 
Constraints on Data Properties 

The goals of all the science programs discussed above 
(and many more, of course) can be accomplished by sat- 
isfying the following minimal constraints^" 

'^^ For a more elaborate listing of various constraints, in- 
cluding detailed specification of various probability distribu- 
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Delivered Image Quality 
(FWHM arcsec) 

Fig. 1. — The image quality distribution measured at the Cerro 
Pachon site using a DIMM (differential image motion monitor) at 
A = 500 nm, and corrected using an outer scale parameter of 30 m 
over an 8.4 m aperture. For details about the outer scale correction 
see Tokovinin (2002). The observed distribution is well described 
by a log-normal distribution, using the parameters shown in the 
figure. 

1. The single visit depth should reach r ^ 24.5. This 
limit is primarily driven by the NEC survey and 
variable sources (e.g., SN, RR Lyrae stars), and 
by proper motion and trigonometric parallax mea- 
surements for stars. Indirectly, it is also driven by 
the requirements on the coadded survey depth and 
the minimum number of exposures placed by WL 
science. 

2. Image quality should maintain the limit set by the 
atmosphere (the median free-air seeing is 0.7 arcsec 
in the r band at the chosen site, see Fig. [T]), and 
not be degraded appreciably by the hardware. In 
addition to stringent constraints from weak lensing, 
good image quality is driven by the required survey 
depth for point sources and by image differencing 
techniques. 

3. Photometric repeatability should achieve 5 mmag 
precision at the bright end, with zeropoint stabil- 
ity across the sky of 10 mmag and band-to-band 
calibration errors not larger than 5 mmag. These 
requirements are driven by the photometric red- 
shift accuracy, the separation of stellar populations, 
detection of low-amplitude variable objects (such 
as eclipsing planetary systems), and the search for 
systematic effects in type la supernova light curves. 

4. Astrometric precision should maintain the limit set 
by the atmosphere, of about 10 mas per visit at the 
bright end (on scales below 20 arcmin). This pre- 
cision is driven by the desire to achieve a proper 
motion accuracy of 0.2 mas/yr and parallax accu- 
racy of 1.0 mas over the course of a 10-year survey 
(see ^2XT|) . 

i :ions, please see the LSST Science R equirements Document 
| |http://www.lsst.org/Science/docs.shtml{ . 
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5. The single visit exposure time (including both ex- 
posures in a visit) should be less than about a 
minute to prevent trailing of fast moving objects 
and to aid control of various systematic effects in- 
duced by the atmosphere. It should be longer than 
~20 seconds to avoid significant efficiency losses 
due to finite readout, slew time, and read noise. 

6. The filter complement should include at least six fil- 
ters in the wavelength range limited by atmospheric 
absorption and silicon detection efficiency (320- 
1050 nm), with roughly rectangular filters and no 
large gaps in the coverage, in order to enable ro- 
bust and accurate photometric redshifts and stellar 
typing. An SDSS-like u band is extremely impor- 
tant for separating low-redshift quasars from hot 
stars, and for estimating the metallicities of F/G 
main sequence stars. A bandpass with an effective 
wavelength of about 1 micron would enable stud- 
ies of sub-stellar objects, high-redshift quasars (to 
redshifts of ~7.5), and regions of the Galaxy that 
are obscured by interstellar dust. 

7. The revisit time distribution should enable determi- 
nation of orbits of Solar System objects and sample 
SN light curves every few days, while accommodat- 
ing constraints set by proper motion and trigono- 
metric parallax measurements. 

8. The total number of visits of any given area of sky, 
when accounting for all filters, should be of the or- 
der of 1,000, as mandated by WL science, the NEO 
survey, and proper motion and trigonometric par- 
allax measurements. Studies of transient sources 
also benefit from a larger number of visits. 

9. The coadded survey depth should reach r ~ 27.5, 
with sufficient signal-to-noise ratio in other bands 
to address both extragalactic and Galactic science 
drivers. 

10. The distribution of visits per filter should enable 
accurate photometric redshifts, separation of stel- 
lar populations, and sufficient depth to enable de- 
tection of faint extremely red sources (e.g., brown 
dwarfs and high-redshift quasars). Detailed sim- 
ulations of photometric redshift estimates sug- 
gest an approximately flat distribution of visits 
among bandpasses (because the system throughput 
and atmospheric properties are wavelength depen- 
dent, the achieved depths are different in different 
bands). The adopted time allocation (see Table 
1) gives a slight preference to the r and i bands 
because of their dominant role in star /galaxy sep- 
aration and weak lensing measurements. 

11. The distribution of visits on the sky should extend 
over at least ^^20,000 deg^ to obtain the required 
number of galaxies for WL studies, with attention 
paid to include "special" regions such as the Eclip- 
tic, Galactic plane, and the Large and Small Mag- 
ellanic Clouds. 

12. Data processing, data products and data access 
should enable efficient science analysis without a 
significant impact on the final uncertainties. To 



TABLE 1 

The LSST Baseline Design and Survey Parameters 



Quantity 


Baseline Design Specification 


Mount Config. 
Final f-Ratio, aperture 
Field of view, etendue 
Plate Scale 
Pixel count 
Wavelength Coverage 
Single visit depths" (5(t) 
Mean number of visits 
Final (coadded) depths" 


"-i— TTi iTTi"!?" rvi rtri 1 fifiri Paul— Ra Vct" 

Alt- azimuth 

f/1.234, 8.4 m 

9.6 deg2, 318 m^deg^ 

50.9 (Um/arcsec (0.2" pix) 

3.2 Gigapix 

320 - 1050 nm, ugrizy 

23.9, 25.0, 24.7, 24.0, 23.3, 22.1 

70, 100, 230, 230, 200, 200 

26.3, 27.5, 27.7, 27.0, 26.2, 24.9 



" The listed values for 5(7 depths in the ugrizy bands, respectively, 
are AB magnitudes, and correspond to point sources and zenith 
observations. See Table 2 for more details. 

enable a fast and efficient response to transient 
sources, the processing latency should be less than 
a minute, with a robust and accurate preliminary 
classification of reported transients. 

Remarkably, even with these joint requirements, mmc 
of the individual science programs is severely over- 
designed. That is, despite their significant scientific di- 
versity, these programs are highly compatible in terms of 
desired data characteristics. Indeed, any one of the four 
main science drivers could be removed, and the remain- 
ing three would still yield very similar requirements for 
most system parameters. As a result, the LSST system 
can adopt a highly efficient survey strategy where a single 
dataset serves most science programs (instead of science- 
specific surveys executed in series). One can think of 
this as massively parallel astrophysics. The vast major- 
ity (about 90%) of the observing time will be devoted to 
a decp-wide-fast survey mode, with the remaining 10% 
of observing time allocated to special programs which 
will also address multiple science goals. Before describ- 
ing these surveys in detail, we discuss the main system 
parameters. 

2.2. The Main System Design Parameters 

Given the minimum science-driven constraints on the 
data properties listed in the previous section, we now 
discuss how they are translated into constraints on the 
main system design parameters: the aperture size, the 
survey lifetime, the optimal exposure time, and the filter 
complement. 

2.2.1. The Aperture Size 

The product of the system's etendue and the survey 
lifetime, for given observing conditions, determines the 
sky area that can be surveyed to a given depth, where the 
etendue is the product of the primary mirror area and the 
field-of-view area. The LSST field-of-view area is maxi- 
mized to its practical limit, 10 deg^, determined by the 
requirement that the delivered image quality be domi- 
nated by atmospheric; seeing at the chosen sitc^ (Ccirro 
Pachon in Northern Chile). A larger field-of-view would 
lead to unacceptable deterioration of the image quality. 
This leaves the primary mirror diameter and survey life- 
time as free parameters. The adopted survey lifetime of 
10 years is a compromise between a shorter time that 
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Fig. 2. — The coadded depth in the r band (AB magnitudes) 
vs. the effective aperture and the survey Hfetime. It is assumed 
that 22% of the total observing time (corrected for weather and 
other losses) is allocated for the r band, and that the ratio of the 
surveyed sky area to the field-of-view area is 2,000. 

leads to an excessively large and expensive mirror (15m 
for a 3 year-long survey and 12m for a 5-year long sur- 
vey), and a smaller telescope that would require more 
time to complete the survey, with the associated increase 
in operations cost. 

The primary mirror size is a function of the required 
survey depth and the desired sky coverage. By and large, 
the anticipated science outcome scales with the num- 
ber of detected sources. For practically all astronomical 
source populations, in order to maximize the number of 
detected sources, it is more advantageous to maximize 
the area first, and then the detection depth"*^ For this 
reason, the sky area for the main survey is maximized 
to its practical limit, 20,000 deg^, determined by the 
requirement to avoid large airmasses {X < 1.5, where 
X = 1/ cos(0) and 9 is the zenith distance), which would 
substantially deteriorate the image quality and the sur- 
vey depth (see eq. ^ . 

With the adopted field-of-view area, the sky coverage 
and the survey lifetime fixed, the primary mirror diam- 
eter is fully driven by the required survey depth. There 
are two depth requirements: the final (coadded) survey 
depth, r ^ 27.5, and the depth of a single visit, r ^ 24.5. 
The two requirements are compatible if the number of 
visits is several hundred (per band), which is in good 
agreement with independent science-driven requirements 
on the latter. 

If the exposure time is doubled and used to double the survey 
area, the number of sources increases by a factor of two. If the 
survey area is kept fixed, the increased exposure time will result 
in ~0.4 mag deeper data (see eq. |6ll. For cumulative source counts 
described by log(Af) = C+k*m, the number of sources will increase 
by more than a factor of two only if fc > 0.75. Apart from low- 
redshift quasars {z < 2), practically all populations have k at most 
0.6 (the so-called Euclidean counts), and faint stars and galaxies 
have k < 0.5. For more details, please see Nemiroff (2003). 



Fig. 3. — The single- visit depth in the r band (5(t detection 
for point sources, AB magnitudes) vs. revisit time, n (days), as a 
function of the effective aperture size. With a coverage of 10,000 
deg^ in two bands, the revisit time directly constrains the visit ex- 
posure time, tyis = 10 n seconds. In addition to direct constraints 
on optimal exposure time, t„is is also driven by requirements on 
the revisit time, n, the total number of visits per sky position over 
the survey lifetime, N^igit, and the survey efficiency, e (see eqs.l- 
3). Note that these constraints result in a fairly narrow range of 
allowed tyis for the main deep-wide-fast survey. 

The required coadded survey depth provides a direct 
constraint, independent of the details of survey execution 
such as the exposure time per visit, on the minimum 
effective primary mirror diameter of 6.5m, as illustrated 
in Fig. H 

2.2.2. The Optimal Exposure Time 

The single visit depth depends on both the primary 
mirror diameter and the chosen exposure time, tvis- In 
turn, the exposure time determines the time interval to 
revisit a given sky position and the total number of visits, 
and each of these quantities has its own science drivers. 
We summarize these simultaneous constraints in terms 
of the single-visit exposure time: 

• The single- visit exposure time should not be longer 
than about a minute to prevent trailing of fast So- 
lar System moving objects, and to enable efficient 
control of atmospheric systematics. 

• The mean revisit time (assuming uniform cadence) 
for a given position on the sky, n, scales as 

days, (1) 



A. 



sky 



10 secy VlO,OOOdeg' 



lOdeg- 

'a 



FOV 



where two visits per night are assumed, and the 
losses for realistic observing conditions have been 
taken into account (with the aid of the Operations 
Simulator described below) . Science drivers such as 
SN and moving objects in the Solar System require 
that n < 4 days, or equivalently tms < 40 seconds 
for the nominal values of Asky and Apov- 



8 



The number of visits to a given position on the sky, 
Nvisit , with losses for reahstic observing conditions 
taken into account, is given by 



/3000\ / T 



\ n ) 



(2) 



The requirement N^sit > 800 again imphes that 
n < 4 and t^^s < 40 seconds if the survey hfetime, 
T 10 years. 

These three requirements place a firm upper limit 
on the optimal visit exposure time of tms < 40 
seconds. Surveying efficiency (the ratio of open- 
shutter time to the total time spent per visit) con- 
siderations place a lower limit on tyis due to fi- 
nite read-out and slew time (the longest acceptable 
read-out time is set to 2 seconds, and the slew and 
settle time is set to 5 seconds, including the read- 
out time for the second exposure in a visit): 



^vis ^~ 9 sec 



(3) 



To maintain efficiency losses below 30% (i.e., at 
least below the limit set by the weather patterns), 
and to minimize the read noise impact, tyis > 20 
seconds. 

Taking these constraints simultaneously into account, 
as summarized in Fig. [3l yielded the following reference 
design: 

1. A primary mirror effective diameter of ^6. 5m. 
With the adopted optical design, described below, 
this effective diameter corresponds to a geometrical 
diameter of ^8m. Motivated by characteristics of 
the existing equipment at the Steward Mirror Lab- 
oratory, which is casting the primary mirror, the 
adopted geometrical diameter is set to 8.4m. 

2. A visit exposure time of 30 seconds (using two 15 
second exposures to efficiently reject cosmic rays), 
yielding e = 77%. 



3. A revisit time of 3 days on average per 10,000 deg^ 
of sky, with two visits per night. 

To summarize, the chosen primary mirror diameter is 
the minimum diameter that simultaneously satisfies the 
depth (r ~ 24.5 for single visit and r ^ 27.5 for coadded 
depth) and cadence (revisit time of 3-4 days, with 30 
seconds per visit) constraints described above. 

2.3. System Design Trade-offs 

We note that the Pan-STARRS project (Kaiser et al. 
2002), with similar science goals as LSST, has adopted 
a distributed aperture design, where the total system 
etendue is a sum of etendue values for an array of small 
telescopes (the prototype PSl telescope has an etendue of 
1/24*'' of the LSST's etendue). Similarly, the LSST sys- 
tem could perhaps be made as two smaller copies with 
6m mirrors, or 4 copies with 4m mirrors, or 16 copies 
with 2m mirrors. Each of these clones would have to 
have its own 3 Gigapixel camera (see below), and given 
the added risk and complexity (e.g., data processing). 



the monolithic design seems advantageous for a system 
with such a large etendue as LSST. 

It is informative to consider the tradeoffs that would 
be required for a system with a smaller aperture, if the 
science requirements were to be maintained. For this 
comparison, we consider a four-telescope version of the 
Pan-STARRS survey (PS4). With an etendue about 6 
times smaller than that of LSST (effective diameters of 
6.5m and 3.0m, and a field-of-view area of 9.6 deg^ vs. 
7.2 deg^), and all observing conditions being equal, the 
PS4 system could in principle use an identical cadence as 
that of LSST. The main difference in the datasets would 
be a faint limit shallower by about 1 mag in a given 
survey lifetime. As a result, for Euclidean populations 
the sample sizes would go down by a factor of 4, while 
for shallower populations (e.g., galaxies around redshift 
of 1) the samples would be smaller by about a factor 2- 
3. The distance limits for nearby sources, such as Milky 
Way stars, would drop to 60% of their corresponding 
LSST values, and the NEO completeness level mandated 
by the U.S. Congress would not be reached. 

If instead the survey coadded depth were to be main- 
tained, then the survey sky area would have to be 6 times 
smaller (^3,500 deg^). If the survey single-visit depth 
were to be maintained, then the exposure time would 
have to be about 6 times longer (ignoring the slight dif- 
ference in the field-of-view area and simply scaling by the 
etendue ratio), resulting in non-negligible trailing losses 
for solar system objects, and either i) a factor of six 
smaller sky area observed within n — Z days, or ii) the 
same sky area revisited every n = 18 days. Given these 
conflicts, one solution would be to split the observing 
time and allocate it to individual specialized programs 
(e.g., large sky area vs. deep coadded data vs. deep 
single-visit data vs. small n data, etc.), as considered by 
the PSl Consortium*^. 

In summary, given the science requirements as stated 
here, there is a minimum etendue of ^300 deg^rr? which 
enables seemingly disparate science goals to he addressed 
with a single data set. A system with a smaller etendue 
would require specialized surveys to address the science 
goals, which results in a loss of surveying efficiency*^. 
The LSST is designed to reach this minimum etendue 
for the science goals stated in its Science Requirements 
Document. 

2.4. The Filter Complement 

The LSST filter complement {ugrizy, see Fig. |4]) is 
modeled after the Sloan Digital Sky Survey (SDSS) sys- 
tem (Fukugita et al. 1996) because of its demonstrated 
success in a wide variety of applications, including pho- 
tometric redshifts of galaxies (Budavari et al. 2003), sep- 
aration of stellar populations (Lenz et al. 1998; Helmi et 
al. 2003), and photometric selection of quasars (Richards 
et al. 2002). The extension of the SDSS system to 
longer wavelengths (the y band at ^^1 micron) is driven 
by the increased effective redshift range achievable with 
the LSST due to deeper imaging, the desire to study 
sub-stellar objects, high-redshift quasars, regions of the 

More information about Pan-STARRS is available from 
http: / /pan-starrs. ifa.hawaii.edu' 

''•^ The converse is also true: for every etendue there is a set of 
optimal science goals that such a system can address with a high 
efficiency. 
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Fig. 4. — The current design of the LSST bandpasses. The 
vertical axis shows the total throughput. The computation includes 
the atmospheric transmission, optics, and the detector sensitivity. 
The detailed design of the y band is still being optimized and the 
figure shows the so-called version which avoids the atmospheric 
water absorption feature at ~0.95 ^m. 

Galaxy that are obscured by interstellar dust, and the 
scientific opportunity enabled by modern CCDs with 
high quantum efhciency in the near infrared. 

The chosen filter complement corresponds to a design 
"sweet spot". We have investigated the possibility of 
replacing the ugrizy system with a filter complement 
that includes only five filters. For example, each filter 
width could be increased by 20% over the same wave- 
length range (neither a shorter wavelength range, nor 
gaps in the wavelength coverage are desirable options), 
but this option is not satisfactory. Placing the red edge 
of the u band blueward of the Balmer break allows op- 
timal separation of stars and quasars, and the telluric 
water absorption feature at 9500 Angstroms effectively 
defines the blue edge of the y band. Of the remaining 
four filters {griz), the g band is already quite wide. As 
a last option, the riz bands could be redesigned as two 
wider bands. However, this option is also undesirable 
because the r and i bands are the primary bands for 
weak lensing studies and for star/galaxy separation, and 
chromatic atmospheric refraction would worsen the point 
spread function for a wider bandpass. 

2.5. The Calibration Methods 

Precise determination of the point spread function 
across each image, accurate photometric and astrometric 
calibration, and continuous monitoring of system perfor- 
mance and observing conditions will be needed to reach 
the full potential of the LSST mission. Extensive pre- 
cursor data including the massive SDSS dataset and our 
own data obtained using telescopes close to the LSST 
site of Cerro Pachon (e.g., the SOAR and Gemini South 
telescopes), as well as telescopes of similar aperture (e.g., 
Subaru), indicate that the photometric and astrometric 
accuracy will be limited not by our instrumentation or 
software, but rather by atmospheric effects. Active op- 
tics will assure superb image quality. 

The required 1% photometric accuracy is driven by our 
requirements on the photometric redshift accuracy, the 
separation of stellar populations, the ability to detect 
low-amplitude variable objects, and the search for sys- 
tematic effects in type la supernova light curves. SDSS 
data taken in good photometric conditions have indeed 
reached the LSST requirements (Padmanabahn et al. 
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Fig. .5. — An example of determination of the atmospheric opac- 
ity by simultaneously fitting a 3-parameter stellar model SED (Ku- 
rucz 1979) and six physical parameters of a sophisticated atmo- 
spheric model (MODTRAN, Anderson et al. 1999) to an observed 
F-type stellar spectrum (Fx). The black line is the observed spec- 
trum and the red line is the best fit. Note that the atmospheric 
water feature around 0.9-1.0 iim is exquisitely well fit. The com- 
ponents of the best-fit atmospheric opacity are shown in Fig. |6] 
Adapted from Burke et al. (2007). 
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Fig. 6. — The components of the best-fit atmospheric opac- 
ity used to model the observed stellar spectrum shown in Fig. [5] 
The atmosphere model (MODTRAN, Anderson et al. 1999) in- 
cludes six components: water vapor (blue), oxygen and other trace 
molecules (green), ozone (red), Rayleigh scattering (cyan), a gray 
term with a transmission of 0.989 (not shown) and an aerosol con- 
tribution proportional to and extinction of 1.3% at A=0.675 
fim (not shown). The black line shows all six components com- 
bined. Adapted from Burke et al. (2007). 



2008), though measurements with ground-based tele- 
scopes typically produce data with errors a factor of two 
or so larger. Analysis of repeated SDSS scans obtained 
in varying observing conditions demonstrates that data 
obtained in traditionally non-photometric conditions can 
also be calibrated with sufficient accuracy (Ivezic et al. 
2007b). 

The LSST calibration plan builds on experience gained 
from the SDSS survey. The process decouples the estab- 
lishment of a stable and uniform internal relative calibra- 
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Fig. 7. — The LSST baseline optical design with its unique 
monolithic mirror: the primary and tertiary mirrors are positioned 
such that they form a continuous compound surface, allowing them 
to be polished into a single substrate. 

tion from the task of assigning absolute optical flux to 
celestial objects. The latter task requires determination 
of a small number of factors appropriate for the accumu- 
lated multi-epoch data set. Celestial sources will be used 
to define the internal photometric system and to monitor 
stability and uniformity of photometric data. There will 
be >100 main-sequence stars with 17 < r < 20 per de- 
tector (14x14 arcmin^) even at high Galactic latitudes. 
Standardization of photometric scales will be achieved 
through direct observation of stars with well-understood 
spectral energy distributions (SEDs). 

While the primary source of data used for photometric 
calibration will be science images taken with the main 
telescope, these images alone will be insufficient to fully 
characterize instrumental and atmospheric temporal and 
spatial variations that can affect photometric measure- 
ments. Auxiliary instrumentation, including a 1.5m cal- 
ibration telescope will provide the calibration parame- 
ters needed for image processing, to calibrate the instru- 
mental response of the LSST hardware (Stubbs & Tonry 
2006), and to measure the atmospheric optical depth as 
a function of wavelength along the LSST line of sight 
(Stubbs et al. 2007; see Figs. [5] and E]). 

2.6. The LSST Reference Design 

We briefly describe the reference design for the main 
LSST system components. Detailed discussion of the 
flow-down from science requirements to system design 
parameters, and extensive system engineering analysis 
can be found in Claver et al. (2008, in prep.). 

2.6.1. Telescope and Site 

The large LSST etendue is achieved in a novel three- 
mirror design (modified Paul-Baker; Angel, Lesser & Sar- 
lot 2000) with a very fast f/1.234 beam. The optical de- 
sign has been optimized to yield a large field of view (9.6 
deg^), with seeing-limited image quality, across a wide 




Fig. 8. — The baseline design (modified three-mirror Paul-Baker) 
for the LSST telescope (current since April 2006). 

wavelength band (320-1050 nm). Incident light is col- 
lected by the primary mirror, which is an annulus with an 
outer diameter of 8.4 m (an effective diameter of 6.5m), 
then reflected to a 3.4m convex secondary, onto a 5m 
concave tertiary, and finally into three refractive lenses 
in a camera (see Fig.[71). This is achieved with an innova- 
tive approach that positions the tertiary mirror inside the 
primary mirror annulus ring, making it possible to fabri- 
cate the mirror pair from a single monolithic blank using 
borosilicate technology. The secondary will be a thin 
meniscus mirror, fabricated from an ultra-low expansion 
material. All three mirrors will be actively supported to 
control wavefront distortions introduced by gravity and 
environmental stresses on the telescope. 

The telescope mount is a compact, stiff structure with 
a fundamental frequency of nearly 10 Hz, which is crucial 
for achieving the required fast slew-and-settle times (see 
Fig. E]). The telescope sits on a concrete pier within a 
carousel dome that is 30 m in diameter (Fig. The 
dome has been designed to reduce dome seeing (local 
air turbulence that can distort images) and to maintain 
a uniform thermal environment over the course of the 
night. The LSST Observatory will be sited atop Cerro 
Pachon in northern Chile, near the Gemini South and 
SOAR telescopes (latitude: S 30° 10' 20.1"; longitude: 
W 70° 48' 0.1"; elevation: 2123 m; the median r band 
zenith seeing: 0.7 arcsec). 

2.6.2. Camera 

The LSST camera provides a 3.2 Gigapixel flat focal 
plane array, tiled by 189 4Kx4K CCD science sensors 
with 10 fim pixels (see Figs. [10] and [TT|). This pixel 
count is a direct consequence of sampling the 9.6 deg^ 
field-of-view (0.64m diameter) with 0.2x0.2 arcsec^ pix- 
els (Nyquist sampling in the best expected seeing of ~0.4 
arcsec). The sensors are deep depleted, back-illuminated 
devices with a highly segmented architecture that enables 
the entire array to be read in 2 seconds. The detectors 
are grouped into 3x3 rafts (see Fig. [12]) , each contain- 
ing its own dedicated front-end and back-end electron- 
ics boards. The rafts are mounted on a silicon carbide 
grid inside a vacuum cryostat, with an intricate thermal 
control system that maintains the CCDs at an operat- 
ing temperature of 180 K. The entrance window to the 
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Fig. 9. — The LSST Observatory: artist's rendering of the dome 
enclosure with the attached summit support building and the LSST 
calibration telescope shown on an adjacent peak (Michael J. Mullen 
Design). 




Fig. 10. — The LSST camera with a person to indicate scale size. 
The camera is positioned in the middle of the telescope and will 
include a filter mechanism and shuttering capability. 
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Fig. 11. — The LSST focal plane. Each cyan square represents 
one 4096 X 4096 pixel sensor. Nine sensors are assembled into a 
raft; the 21 rafts are outlined in red. There are 189 science sensors, 
each with 16.8 Mpix, for a total pixel count of 3.2 Gpix. 



cryostat is the third of the three refractive lenses in the 
camera. The other two lenses are mounted in an optics 
structure at the front of the camera body, which also 
contains a mechanical shutter, and a carousel assembly 



Imaging Sensors 




Fig. 12. — The LSST raft module with integrated front-end 
electronics and thermal connections. Each raft includes 9 sensors 
and can be replaced on a short timescale. 
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Fig. 13. — The three-layered architecture of the data manage- 
ment system (application, middleware and infrastructure layers) 
enables scalability, reliability, and evolutionary capability. 

that holds five large optical filters. A sixth optical filter 
will also be fabricated and can replace any of the five via 
a procedure accomplished during daylight hours. 



2.6.3. Data Management 

The rapid cadence of the LSST observing program will 
produce an enormous volume of data, ~30 TB per night, 
leading to a total database over the ten years of oper- 
ations of 60 PB for the raw data, and 30 PB for the 
catalog database. The total data volume after process- 
ing will be several hundred PB, processed using ^100 
TFlops of computing power. Processing such a large vol- 
ume of data, converting the raw images into a faithful 
representation of the Universe, automated data quality 
assessment, and archiving the results in useful form for 
a broad community of users is a major challenge. 

The data management system is configured in three 
levels: an infrastructure layer consisting of the comput- 
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Fig. 14. — The application layer of the data management system. 
The arrows indicate direction of the data flow. The main pipelines 
enable nightly processing, calibration, and data releases 

ing, storage, and networking hardware and system soft- 
ware; a middleware layer, which handles distributed pro- 
cessing, data access, user interface, and system opera- 
tions services; and an applications layer, which includes 
the data pipelines and products and the science data 
archives (see Fig. lT^ . The application layer (see Fig. rH)) 
is organized around the data products being produced. 
The nightly pipelines (see Fig. [T5|) are based on image 
subtraction, and are designed to rapidly detect interest- 
ing transient events in the image stream and send out 
alerts to the community within 30 seconds from com- 
pleting the image readout. The data release pipelines 
(see Fig. ITB)) . in contrast, arc intended to produce the 
most completely analyzed data products of the survey, 
in particular those that measure very faint objects and 
cover long time scales. A new run begins each year, pro- 
cessing the entire survey data set that is available. The 
data release pipelines consume most of the computing 
power of the data management system. The calibration 
products pipeline produces the wide variety of calibra- 
tion data required by the other pipelines. All of these 
pipelines are architected to make efficient use of linux 
clusters with thousands of nodes. 

There will be both mountain summit and base com- 
puting facilities, as well as a central archive facility and 
multiple data access centers. The data will be trans- 
ported over existing high-speed optical fiber links from 
South America to the U.S. (see Fig. [TT]). Although the 
data processing center will have substantial computing 
power (~100 TFlops, equal to the world's most power- 
ful computer in 2004), the continuation of current trends 
suggests that it will not even qualify for the top 500 list 
by the time of first light in 2014. Hence, while LSST 
is making a novel use of advances in information tech- 
nology, it is not taking the risk of pushing the expected 
technology to the limit. 

3. ANTICIPATED DATA PRODUCTS AND THEIR 
CHARACTERISTICS 

The LSST observing strategy is designed to maximize 
the scientific throughput by minimizing slew and other 
downtime and by making appropriate choices of the filter 
bands given the real-time weather conditions. We first 
describe a simulator that has been developed to evaluate 
this process and then illustrate with two examples how 
its outputs have been used to predict LSST performance. 

3.1. The LSST Operations Simulator 

The LSST Operations Simulator was developed to en- 
able a detailed quantitative analysis of the various science 
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Fig. 15. — The structure of nightly pipchncs. Their main task 
is to compare new data to the previously accumulated data and 
search for transient, variable and moving objects. 




Fig. 16. — The structure of data release pipelines. These 
pipelines are intended to produce the most completely analyzed 
data products of the survey, and will consume most of the com- 
puting power of the data management system. 
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tradeoffs described in this paper. It contains detailed 
models of site conditions, hardware and software per- 
formance, and an algorithm for scheduling observations 
which will, eventually, drive the largely robotic observa- 
tory. Observing conditions include a model for seeing 
derived from an extensive body of on-site MASS/DIMM 
(Multi- Aperture Scintillation Sensor and Differential Im- 
age Motion Monitor) measurements obtained during site 
selection and characterization (see Fig. [J). It not only 
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reproduces the observed seeing distribution, but includes 
the auto-correlation spectrum of seeing with time over in- 
tervals from minutes to seasons. Weather data are taken 
from ten years of hourly measurements at nearby Cerro 
Tololo. The time history of site conditions is important 
if the simulation is to faithfully reproduce a sequence 
of observations. Thus, for example, the simulator cor- 
rectly represents the variation of limiting magnitude be- 
tween pairs of observations used to detect NEOs and the 
correlation between, for example, seasonal weather pat- 
terns and observing conditions at any given point on the 
sky. The signal-to-noise ratio of each observation is de- 
termined using a sky background model which includes 
the dark sky brightness in each filter, the effects of seeing 
and atmospheric transparency, and a detailed model for 
scattered light from the moon and/or twilight at each 
observation. The time taken to move from one obser- 
vation to the next is given by a detailed model of the 
camera, telescope, and dome. It includes such effects as 
the acceleration/deceleration profiles employed in mov- 
ing in altitude, azimuth, camera rotator, dome azimuth, 
and wind/stray light screen altitude, the time taken to 
damp vibrations excited by each slew, cable wrap, time 
taken for active optics lock and correction as a function 
of slew distance, filter change, and focal plane readout. 

Observations are scheduled by a ranking algorithm. 
After a given exposure, all possible next observations 
are assigned a score which depends upon their locations, 
times, and filters according to a set of scientific require- 
ments which can vary with time and location. For exam- 
ple, if an ecliptic field has been observed in the r band, 
the score for another r-band observation of the same field 
will initially be quite low, but it will rise in time to peak 
about an hour after the first observation, and decline 
thereafter. This results in these observations being ac- 
quired as pairs roughly an hour apart, which enables ef- 
ficient association of NEO detections. To ensure uniform 
sky coverage, fields with fewer previous observations will 
be scored more highly than those which have already 
been observed more frequently. 

Once all possible next observations have been scored 
for scientific priority, their scores are modified according 
to observing conditions (e.g., seeing, airmass, and sky 
brightness) and to criteria such as low slew time to move 
from the current position, time spent to change filters, 
etc. The highest-ranked observation is then performed, 
and the cycle repeats. The result of a simulator run 
is a detailed history of which locations on the sky were 
observed when, in what filter, and with what sky back- 
ground, seeing and other observing conditions. It takes 
about a week to produce a decade-long simulation using 
an average PC. 

3.2. The Baseline LSST Surveys 

The fundamental basis of the LSST concept is to scan 
the sky deep, wide, and fast, and to obtain a dataset 
that simultaneously satisfies the majority of the science 
goals. This concept, the so-called "universal cadence", 
will yield the main deep-wide-fast survey (typical single 
visit depth of r ^ 24.5) and use about 90% of the ob- 
serving time. The remaining 10% of observing time will 
be used to obtain improved coverage of parameter space 
such as very deep (r ~ 26) observations, observations 
with very short revisit times (~1 minute), and observa- 
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Fig. 18. — The distribution of the r band visits on the sky 
for one simulated reaUzation of the baseUne main survey. The 
sky is shown in Aitoff projection in equatorial coordinates and the 
number of visits for a 10-year survey is color-coded according to the 
inset. The two regions with smaller number of visits than the main 
survey ("mini-surveys") are the Galactic plane (arc on the left) 
and the so-called "northern Ecliptic region" (upper right). The 
region around the South Celestial Pole will also receive substantial 
coverage (not shown here). 

tions of "special" regions such as the Ecliptic, Galactic 
plane, and the Large and Small Magellanic Clouds. We 
are also considering a third type of survey, micro-surveys, 
that would use about 1% of the time (which still repre- 
sents 25 nights on a unique 8m-class telescope). 

3.2.1. The Main Deep-Wide-Fast Survey 

The observing strategy for the main survey will be opti- 
mized for the homogeneity of depth and number of visits. 
In times of good seeing and at low airmass, preference is 
given to r-band and i-band observations. As often as 
possible, each field will be observed twice, with visits 
separated by 15-60 minutes. This strategy will provide 
motion vectors to link detections of moving objects in 
the Solar System, and fine-time sampling for measuring 
short-period variability. The ranking criteria also ensure 
that the visits to each field are widely distributed in po- 
sition angle on the sky and rotation angle of the camera 
in order to minimize systematic effects in galaxy shape 
determination. 

The universal cadence provides most of LSST's power 
for detecting Near Earth Objects (NEO) and naturally 
incorporates the southern half of the ecliptic within its 
20,000 square degrees (the northern half lies above the 
desired airmass limits, X < 1.5). NEO survey complete- 
ness for the smallest bodies (~140m in diameter, per the 
Congressional NEO mandate) is greatly enhanced, how- 
ever, by the addition of a crescent on the sky within 10 
degrees of the northern ecliptic (see Fig. [T8|) . Thus, the 
"northern Ecliptic proposal" extends the universal ca- 
dence to this region using the r and i filters only, along 
with more relaxed limits on airmass and seeing. Relaxed 
limits on airmass and seeing are also adopted for '^700 
deg^ around the south celestial pole, allowing coverage 
of the Large and Small Magellanic Clouds. 

Finally, the universal cadence proposal excludes obser- 
vations in a region of 1,000 square degrees around the 
Galactic Center, where the high stellar density leads to a 
confusion limit at much brighter magnitudes than those 
attained in the rest of the survey. Within this region, 
the Galactic Center proposal provides 30 observations in 
each of the six filters, distributed roughly logarithmically 
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in time (it may not be necessary to use the bluest u and 
g filters for this heavily extincted region). 

The resulting sky coverage for the LSST baseline ca- 
dence, based on detailed operations simulations, is shown 
for the r band in Fig. 1181 The anticipated total number 
of visits for a ten-year LSST survey is about 2.8 million 
(~5.6 million 15-second long exposures). The per-band 
allocation of these visits is shown in Table 1. 

3.2.2. Mini-surveys 

Although the uniform treatment of the sky provided by 
the universal cadence proposal can satisfy the majority 
of LSST scientific goals, roughly 10% of the time will be 
allocated to other strategies that significantly enhance 
the scientific return. These surveys aim to extend the 
parameter space accessible to the main survey by going 
deeper or by employing different time/filter sampling. 

As an example of a mini-survey, consider a program 
that uses one hour of observing time per night to observe 
a relatively small region of sky to substantially greater 
depth in individual visits. Accounting for read-out time 
and filter changes, it could obtain about 50 consecutive 
15-second exposures in each of four filters in an hour. 
If a field is visited every two days over four months, 
about 600 deg^ can be observed with this cadence over 
10 years. Taking weather into account, the selected fields 
would each have on average about 40 hour-long sequences 
of 200 exposures each. Each observation in a sequence 
would have an equivalent 5-cr depth of r ^--^ 24, and each 
filter subsequence when coadded would be 2 magnitudes 
deeper than the main survey visits (r ~ 26.5). When 
all 40 sequences and the main survey visits are coadded, 
they would extend the depth to r ^ 28. 

This data set would be excellent for a wide variety of 
science programs. The individual sequences would be 
sensitive to 1% variability on hourly time scales, allow- 
ing discovery of planetary eclipses. If these fields were 
selected at Galactic latitudes of |6| ~ 30 deg, they would 
include about 10 million stars observed with signal-to- 
noise ratio above 100 in each visit (stars brighter than 
r ~ 21). When subsequences from a given night were 
coadded, they would provide dense time sampling to a 
faint limit (in 4 bands, every 2 days over 120 days, mod- 
ulo weather, to r 26.5), and thus enable deep searches 
for SN, trans-Neptunian objects, and other faint tran- 
sient, moving and variable sources. For example, the SN 
sample would be extended to redshifts oi z ^ 1.2, with 
more densely sampled light curves than obtained from 
the universal cadence. Such sequences would also serve 
as excellent tests of our calibration procedures. 

The baseline universal cadence is by no means the 
definitive plan for the entire survey. Rather, it repre- 
sents a proof of concept that it is indeed possible to 
design a universal cadence which addresses a wide va- 
riety of science goals in a nearly optimal way. We are 
undertaking a vigorous and systematic research effort to 
explore the enormously large parameter space of possible 
surveys. The commissioning period will be used to test 
the usefulness of various observing modes and to explore 
alternative strategies. Proposals from the community, 
through the science collaborations (see § [5]) , for special- 
ized cadences (such as mini-surveys and micro-surveys) 
will also be considered. 



TABLE 2 
The Parameters From Eqs. [5]and[6] 
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21.8 


22.0 


21.3 
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19.1 
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0.77 


0.73 


0.70 


0.67 


0.65 
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0.038 


0.039 


0.039 


0.040 


0.040 




23.60 


24.57 


24.57 


24.47 


24.19 


23.74 




0.48 


0.21 


0.10 


0.07 


0.06 


0.06 
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23.9 


25.0 


24.7 


24.0 


23.3 


22.1 




0.21 


0.16 


0.14 


0.13 


0.13 


0.13 



" The expected median zenith sky brightness at Cerro Pachon, 
assuming mean solar cycle and 3-day old moon (mag/arcsec^). 
^ The expected delivered median zenith seeing (arcsec). For larger 
airmass, X, seeing is proportional to X^''-^ . 

The band-dependent parameter from Eq. [5] 

The band-dependent parameter from Eq. [6] 

Adopted atmospheric extinction. 
^ The typical 5(t depth for point sources at zenith, assuming 
exposure time of 2x15 sec, and observing conditions as listed. For 
larger airmass the 5(t depth is brighter; see the bottom row. 
^ The loss of depth at the median airmass oi X = 1.2 due to 
seeing degradation and increased atmospheric extinction. 

3.3. Detailed Analysis of Simulated Surveys 

As examples of analysis enabled by the Operations 
Simulator, we describe determination of the complete- 
ness of the LSST NEO survey, and estimation of er- 
rors expected for trigonometric parallax and proper mo- 
tion measurements. In both examples, the conclusions 
crucially depend on assumed signal-to-noise ratios, de- 
scribed next. 

3.3.1. Expected Signal-to-Noise Ratio 

The output of operations simulations is a data stream 
consisting of a position on the sky and the time of obser- 
vation, together with observing conditions such as seeing 
and sky brightness. The expected photometric error (the 
inverse of signal-to-noise ratio) for a single visit can be 
written as 

C^l = f^sys + errand' (4) 

where Urand is the random photometric error and Ugys 
is the systematic photometric error (includes errors due 
to, e.g., imperfect modeling of the point spread func- 
tion, but does not include errors in absolute photometric 
zeropoint). The calibration system and procedures are 
designed to maintain agys < 0.005 mag. Based on SDSS 
experience (Sesar et al. 2007), the random photometric 
error, as a function of magnitude, is well described^^ by 

^rand = (0.04-7)x + 7x2 (mag2), (5) 

with X = 100-4 (m-ms)^ Here mg is the 5a depth (for 
point sources) in a given band, and 7 depends on the sky 
brightness, readout noise, etc. Using the LSST exposure 
time calculator's (Gee et al. 2007), we have obtained 
the values of 7 listed in Table 2. The 5cr depth for point 
sources is determined from 

TO5 = C,n + 0.50 {m,ky - 21) + 2.5 logio(O.7/0) + 

+1.25 1ogio(i„,/30)-fc„(X-l) (6) 

Eq.[5]can be derived from (T^^nd = tV/S, where N is noise and 
S is signal, and by assuming that A^^ = + aS. The constants 
No and a can be expressed in terms of a single unknown constant 
7 by using the co ndition that errand = 0.2 for m = m s. 
*^ Available at |http: / /dls. physics. ucdavis.edu/etc| 
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where nisky is the sky brightness (mag/arcsec^), 6 is the 
seeing (FWHM, in arcsec), tvis is the exposure time (sec- 
onds), k is the atmospheric extinction coefficient, and 
X is airmass. The constants Cm depend on the over- 
all throughput of the instrument and are determined us- 
ing the LSST exposure time calculator. The assump- 
tions built into the calculator were tested using SDSS 
observations^^, and by comparing the predicted depths 
to the published performance of the Subaru telescope 
(Kashikawa et al. 2003). The adopted values for Cm 
and k are listed in Table 2, as well as expected ms in 
nominal observing conditions. 

3.3.2. The NEO Completeness Analysis 

To assess the LSST completeness for NEOs (more pre- 
cisely, for a 20% subset of NEOs, the potentially haz- 
ardous asteroids or PHAs), the PHA population is rep- 
resented by a size-limited complete sample of 800 true 
PHAs whose orbital elements are taken from the Minor 
Planet Center. The simulated baseline survey is used 
to determine which PHAs are present in each exposure 
and at what signal-to-noise ratio they were observed. In 
addition to seeing, atmospheric transparency, and sky 
background effects (see eq. IH]), the signal-to- noise com- 
putation takes into account losses due to non-optimal 
filters and object trailing. Using SDSS observations of 
asteroids (Ivezic et al. 2001), we adopt the following 
mean colors to transform limiting (AB) magnitudes in 
LSST bandpasses (listed in Table 2) to an 'effective" 
limiting magnitude in the standard V band: V ~ m = 
(-2.1, -0.5, 0.2, 0.4, 0.6, 0.6) for m = (u, g, r, i, z, y). Due 
to very red V — u colors, and the relatively bright lim- 
iting magnitude in the y band, the smallest objects are 
preferentially detected in the griz bands. The correc- 
tion for trailing is implemented by subtracting from the 
right-hand side of cq. [6] 



Am 



traili] 



1.25 logio(l 4-0.0267- 



(7) 



where the object's velocity, v, is expressed in deg/day. 
For the nominal exposure time of 30 seconds and 9 — 0.7 
arcsec, the loss of limiting magnitude is 0.16 mag for 
V = 0.25 deg/day, typical for objects in the main asteroid 
belt, and 0.50 mag for v = 1.0 deg/day, typical of NEOs 
passing near Earth. 

Above 60% completeness, one magnitude of depth is 
roughly equivalent to 10% in completeness. An object's 
orbit is considered to be cataloged if the object was de- 
tected on at least three nights during a single lunation, 
with a minimum of two visits per night. The same cri- 
terion was used in NASA studies''^, and is confirmed as 
reliable by a detailed analysis of orbital linking and de- 
termination using the MOPS code (Jedicke et al. 2005) 
developed by the Pan-STARRS project (and adopted by 

*® SDSS data typically reach a 5cr depth for point sources of 
r = 22.5 with an effective aperture of D = 2.22m, an exposure time 
of tiiis = 54 sec, the median r band sky brightness of r^^y = 20.9 
mag/arcsec-^, the median seeing o{ 6 = 1.5 arcsec, and the median 
airmass of X = 1.3. In comparison, the LSST loss of depth is 0.32 
mag due to shorter exposures, while the gains are 1.17 mag due 
to larger aperture, 0.83 mag due to better seeing, 0.20 mag due 
to fainter sky, and 0.10 mag due to better throughput, for the net 
gain of ^^2.0 mag. 

The NASA 2007 NEO study is available from 
|http://neo.jpl.nasa.gov/neo/report2007.html| 
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Fig. 19. — Completeness of the LSST survey for PHAs brighter 
than a given absolute magnitude (related to the size of the object 
and albedo; H=22 mag is equivalent to a typical 140m asteroid 
and H=2A mag is equivalent to a 50m asteroid). Two scenarios 
are shown: the lower curve is the 10-year long baseline survey 
with a 5% NEO optimization, and it reaches a completeness of 
84%. The upper dashed curve results from spending 15% of the 
observing time in an NEO optimized mode, and running the survey 
for 12 years. It meets the 90% completeness level for 140m objects 
mandated by the Congress. 

LSST in a collaborative effort with Pan-STARRS). The 
MOPS software system and its algorithms are signifi- 
cantly more advanced than anything fielded for this pur- 
pose to date. Realistic MOPS simulations show >99% 
linking efficiency across all classes of solar system ob- 
jects. For the LSST baseline cadence, objects counted as 
cataloged are observed on 20 different nights on average 
over ten years. A more stringent requirement that an 
object must be detected on at least five nights decreases 
the completeness by typically 3%. The completeness is 
also a function of the assumed size distribution of NEOs: 
the flatter the distribution, the higher the completeness. 
If the latest results for the NEO size distribution by A. 
Harris (priv. comm.) are taken into account, the com- 
pleteness increases by 1-2%. Due to these issues, the 
completeness estimates have a systematic uncertainty of 
at least 2%. 

The LSST baseline cadence provides orbits for 82% of 
PHAs larger than 140 meters after 10 years of operations. 
With a minor change of this cadence, such as requiring 
that all observations in the so-called NE region (North 
Ecliptic region, defined by S > 5) are obtained in the 
r band, the completeness for 140m and larger PHAs is 
84%, with 90% completeness reached for 200m and larger 
objects. The completeness curve as a function of an ob- 
ject's size is shown in Fig. [19] (lower curve). This cadence 
spends 5% of the total observing time on NEO-optimized 
observations in the NE region. 

Various adjustments to the baseline cadence can boost 
the completeness for 140m and larger PHAs to 90%. We 
find that such variations can have an unacceptably large 
impact on other science programs, if the 90% complete- 
ness is to be reached within the 10 year survey lifetime. 
However, with a minor adjustment of the baseline ca- 
dence, such that 15% of the time is spent in the NE re- 
gion to reach fainter limiting magnitudes, this complete- 
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ness level can be reached with a 12-year long survey, and 
with a negligible effect on other science goals. The com- 
pleteness curve as a function of an object's size for such 
a modified cadence is shown in Fig. [19] (upper curve). 

Our analysis assumes that no NEOs are known prior 
to LSST. Currently known NEOs do not have a signifi- 
cant impact on this calculation. However, if a precursor 
survey, such as Pan-STARRS 4, operated for three years 
prior to LSST, the time to fulfill the Congressional man- 
date by LSST could be shortened by about a year (the 
current completeness level does not have a significant im- 
pact). 

3.3.3. The Expected Accuracy of Trigonometric Parallax 
and Proper Motion Measurements 

Given the observing sequence for each sky position in 
the main survey, we generate a time sequence of mock as- 
trometric measurements. The assumed astrometric accu- 
racy is a function of signal-to-noise ratio, SNR. Random 
astrometric errors per visit are modeled as 9/ SNR, with 
9 ~ 700 mas and SNR determined using eq. \6\ The esti- 
mated proper motion and parallax accuracy at the bright 
end (r < 20) is driven by systematic errors due to the 
atmosphere. Systematic errors of 10 mas are added in 
quadrature, and are assumed to be uncorrelated between 
different observations of a given object. Systematic and 
random errors become similar at about r — 22, and there 
are about 100 stars per LSST sensor (0.05 sq. deg) to 
this depth (and fainter than the LSST saturation limit 
at r ~ 16) even at the Galactic poles. 

Precursor data from the Subaru telescope indicate that 
systematic errors of 10 mas on spatial scales of several 
arcmins are realistic. Even a drift-scanning survey such 
as SDSS delivers uncorrelated systematic errors (domi- 
nated by seeing effects) at the level of 20-30 mas (mea- 
sured from repeated scans, Pier et al. 2003), and the 
expected image quality for LSST will be twice as good 
as for SDSS. Furthermore, there are close to 1000 galax- 
ies per sensor with r < 22, which will provide exquisite 
control of systematic astrometric errors as a function of 
magnitude, color and other parameters, and thus enable 
absolute proper motion measurements. 

The astrometric transformations for a given CCD and 
exposure, and proper motion and parallax for all the 
stars from a given CCD, are simultaneously solved for 
using an iterative algorithm. The astrometric transfor- 
mations from pixel to sky coordinates are modeled using 
low-order polynomials and standard techniques devel- 
oped at the U.S. Naval Observatory (Monet et al. 2003). 
The expected proper motion and parallax errors for a 
10-year long baseline survey, as a function of apparent 
magnitude, are summarized in Table 3. Blue stars (e.g., 
F/G stars) fainter than r 23 will have about 50% larger 
proper motion and parallax errors than given in the table 
due to decreased numbers of z and y band detections. 
The impact on red stars is smaller due to a relatively 
small number of observations in the u and g bands, but 
extremely red objects, such as L and T dwarfs, will def- 
initely have larger errors, depending on details of their 
spectral energy distribution. After the first three years of 
the survey, the proper motion errors are about five times 
as large, and parallax errors are about twice as large as 
the values given in Table 3; the errors scale as t~^^^ and 
^-1/2^ respectively. This error behavior is a strong in- 



TABLE 3 

The expected proper motion, parallax and accuracy for a 
10-year long baseline survey. 
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mag 


mas 


mas 


mas / yr 


4 
mag 


mag 


21 


11 


0.6 


0.2 


0.01 


0.005 


22 


15 


0.8 


0.3 


0.02 


0.005 


23 


31 


1.3 


0.5 


0.04 


0.006 


24 


74 


2.9 


1.0 


0.10 


0.009 



" Typical astrometric accuracy (rms per coordinate per visit); 
' Parallax accuracy for 10-year long survey; 

Proper motion accuracy for 10-year long survey; 

Photometric error for a single visit (two 15-second exposures); 

Photometric error for coadded observations (see Table 1). 

dependent argument for a survey lifetime of at least 10 
years (c.f. §2). 

For comparison with Table 3, the SDSS-POSS proper 
motion measurements have an accuracy of ~5 mas/yr 
per coordinate at r = 20 (Munn et al. 2004). Gaia is 
expected to deliver parallax errors of 0.3 mas and proper 
motion errors of 0.2 mas/yr at its faint end at r ~ 20. 
Hence, LSST will smoothly extend Gala's error vs. mag- 
nitude curve 4 magnitudes fainter. 

3.4. Data Products 
3.4.1. Standard processing and data products 

The LSST data system is being designed to enable as 
wide a range of science as possible. Standard data prod- 
ucts, including calibrated images and catalogs of detected 
objects and their attributes, will be provided both for in- 
dividual exposures and the deep incremental data coad- 
dition. About 2 billion objects will be routinely moni- 
tored for photometric and astrometric changes, and any 
transient events (non-recurrent objects with statistically 
significant photometric change; about 10,000 per night 
on average) will be posted in less than 60 seconds (with a 
goal of 30 seconds) via web portals including the Virtual 
Observatory. For the "static" sky, there will be yearly 
database releases listing many attributes for billions of 
objects. This database will grow in size to about 30 PB 
and about 20 billion objects^® in ten years, and will in- 
clude other metadata (parameter error estimates, system 
data, seeing summary, etc). For a comparison, the SDSS 
Data Release 6 includes about 287 million unique objects 
(Adelman-McCarthy et al. 2006). 

The SDSS experience is that the distribution of sizes of 
data portal requests follows a power law, ranging from 
frequent simple queries or image cutout downloads to 
infrequent major downloads or large calculations. The 
same distribution is expected for LSST, and we are de- 
signing the system to enable this entire distribution. This 
data archive server will involve different resources de- 
ployed in different ways. We plan to optimize the mix by 
shipping algorithms to the data for those cases that need 
it (and offering local parallel computation) and transfer- 
ring the data to the user for other use cases. While the 
data archive will be at the National Center for Super- 
computing Applications (University of Illinois), we plan 

It is assumed that about 10 billion stars and 10 billion galax- 
ies will be detected by LSST. The estimate of about 10 billion 
detected stars is sensitive to assumptions about the distribution 
of the interstellar dust in the Galactic plane, and cadence details. 
Plausible estimates range from 5 billion to 20 billion stars. 
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multiple data portals for community access to both the 
pixel data and the database. We are designing for a wide 
range of queries and access patterns in order to serve the 
broadest range of scientists as well as interested ama- 
teurs. The LSST database will be interoperable with 
those of surveys in other wavebands, via the Virtual Ob- 
servatory. 

3.4.2. Non-standard data processing and data access 

The database is designed to allow evolution in the na- 
ture of the science questions that are asked of it. This 
design is more than simply configuring it to allow for ex- 
pected growth in query complexity with time. For exam- 
ple, we anticipate that some users will submit queries on 
the images themselves. We plan for two pathways, out- 
side of the database query system, for non-standard pixel 
processing. The first is through a pipeline interface that 
would allow a user group to attach their own pipeline 
processing system at a Data Access Center. This ap- 
proach is appropriate for major undertakings such as an 
independent weak lensing analysis starting at the image 
level, or new ways to coadd the data. The second ap- 
proach is simply to include the desired processing, such 
as a new quantity to be included in the Object database, 
within the data release processing, and pick it up on the 
next data release run. Much of the science will be statis- 
tical in nature and will involve large correlation calcula- 
tions with user-specified precision. We plan on enabling 
open source community driven applications development. 
Finally, special data pipelines developed by the Science 
Collaborations (see § [5]) will populate new parameters in 
the database. 

4. EXAMPLES OF LSST SCIENCE PROJECTS 

The design and optimization of the LSST system lever- 
ages its unique capability to scan a large sky area to a 
faint flux limit in a short amount of time. The main 
product of the LSST system will be a multi-color ugrizy 
image of about half the sky to unprecedented depth 
(r ^ 27.5). For a comparison, the best analogous con- 
temporary dataset is that of SDSS, which provides ugriz 
images of about a quarter of the sky to r ~ 22.5, with 
twice as large seeing (see Figs. [20l and [21]) . A major ad- 
vantage of LSST is the fact that this deep sky map will 
be produced by taking hundreds of shorter exposures (see 
Table 1). Each sky position within the main survey area 
will be observed about 1000 times, with time scales span- 
ning seven orders of magnitude (from 30 sec to 10 years), 
and produce over a trillion photometric measures of ce- 
lestial sources. 

It is not possible to predict all the science that LSST 
data will enable. We now briefly discuss a few projects 
to give a flavor of anticipated studies, and organize them 
by the four science themes that drive the LSST design 
(although some projects span more than one theme). 

4.1. Probing Dark Energy and Dark Matter 

A unique aspect of LSST as a probe of dark energy 
and matter is the use of multiple cross-checking probes 
that reach unprecedented precision (see Fig. [22]) . The 
joint analysis of LSST weak lensing and BAO is partic- 
ularly powerful in constraining the dynamical behavior 




Fig. 20. — A comparison of ~ 7.5 x 7.5 arcmin^ images of the 
same area of sky (centered on q=9'' 20' 47" and <5=30° 8' 12") 
obtained by the SDSS (top, r < 22.5) and the Deep Lens Survey 
(bottom, r < 24.5). The depth gain for the bottom image is mostly 
due to the lower surface brightness limit, which is also responsible 
for the apparent increase of galaxy sizes. LSST will obtain ~100 
gri color images (~200 for the riz composites) of each point over 
half the Celestial sphere (20,000 deg'^, equivalent to 1.3 million 
^ 7.5 X 7.5 arcmin^ regions), and with better seeing. After their 
coaddition, the final image will be another ~ 3 mag deeper (a faint 
limit of r = 27.5 for point sources). 

of dark energy, i.e., how it evolves with cosmic time or 
redshift (Hu & Jain 2004; Zhan 2006). By simultane- 
ously measuring the growth of large-scale structure, and 
luminosity and angular distances as functions of redshift 
(via weak lensing, BAO, SN, and cluster counting), LSST 
data will reveal whether the recent cosmic acceleration 
is due to dark energy or modified gravity (Lue, Scoc- 
cimarro & Starkman 2004; Knox, Song & Tyson 2006; 
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Fig. 21. — A comparison of angular resolution for images ob- 
tained by the SDSS (top, median seeing of 1.5 arcsec) and expected 
from LSST (bottom, seeing of 0.7 arcsec). The images show a 
lensed SDSS quasar (SDSS J1332+0347, Morokuma et al. 2007), 
and the bottom image was taken with Suprime-cam at Subaru. 
Adapted from Blandford et al. (2007). 

Ishak, Upadhye & Spergel 2006; Jain & Zhang 2007). 
Over a broad range of accessible redshifts, the simple lin- 
ear model for the dark energy equation of state is a poor 
representation of more general dark energy theories. Al- 
brecht and Bernstein (2007) and Barnard etal. (2008) 
have recently shown that in a high-dimensional dark en- 
ergy model space, LSST data could lead to a hundred- 
to thousand-fold increase in precision over precursor ex- 
periments about to be undertaken, thereby confirming 
its status as a premier Stage IV experiments^. 

The power and accuracy of LSST dark energy and 
dark matter probes is derived from samples of several 
billion galaxies and millions of Type la supernovae. The 

See the Dark Energy Task Force report (Albrecht et al. 2006) 
for a definition of Stage IV experiments. 
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Fig. 22. — The top panel shows forecasts of LSST errors on 
the dark energy equation-of-state parameters wq and Wa {w{z) = 
P/p = uiQ -I- (1 — a)wa, where a = (I + for measurements 

based on baryon acoustic oscillations, weak lensing, supernovae 
and the three combined. The constraints are marginalized over 
9 other cosmological parameters including the curvature and over 
120 parameters that model the linear galaxy clustering bias, pho- 
tometric redshift bias, and rms photometric redshift error (Zhan 
2006). Simple additive systematic errors of the shear and galaxy 
power spectra are also included. The joint analysis of BAO and 
WL significantly reduces the error ellipse area because the error 
distribution of photometric redshifts can be calibrated with galaxy 
power spectra, and the galaxy bias can be determined with galaxy 
and WL shear power spectra. The bottom panel is analogous to 
the top panel, but for the physical matter density ujm = ^m.h'^ and 
curvature term $7^,. Combining all of LSST's cosmological probes 
further increases the precision. Different theoretical models of dark 
energy span this WQ-Wa space. 

nominal high-SNR sample defined by z < 25 (SNR>25 
for point sources) will include four billion galaxies (56 
arcmin^^) with a mean photometric redshift accuracy 
of 1-2% (relative error for 1 -I- z), and with only 10% 
of the sample with errors larger than 4%. The median 
redshift for this sample will be z ~1.2, with the third 
quartile at z ~ 2. For a subsample of 2 billion galaxies 
further constrained by flux limits in the g and z bands, 
the photometric redshift errors will be two to three times 
smaller. It may be possible to further improve photomet- 
ric redshift calibration using non-conventional methods 
(Newman 2008). The surface density of galaxies useful 
for weak lensing analysis will be about 40 arcmin"^. 
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Fig. 23. — The lensing shear auto power spectra constructed 
from 5 redshift bins (Z is the multipole moment of the distribution 
on the sky). Only the 5 auto-power spectra of each redshift bin 
among the available 15 cospectra are displayed, and the solid curves 
show the predictions for the concordance ACDM model. The boxes 
show the expected l-cr measurement error {Al/l = 0.2) from the 
full LSST 10-year survey due to the sample (i.e., cosmic) variance 
which dominates at about I < 1000, and intrinsic ellipticities which 
dominate at I > 1000 (Zhan 2006). 
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Fig. 25. — Marginalized l-cr errors on the comoving distances 
from LSST WL (solid line and open circles), BAO (dashed line and 
open squares), and joint WL and BAO (dotted line). The distance 
parameters, Di. . . D14, are evenly spaced in log(l+2) from 2i=0.14 
to 214=5. Lines represent results with cosmological parameters 
floating, i.e., marginalized with weaker-than-5-year-WMAP priors, 
and symbols are for results with Planck priors (Zhan, Knox & 
Tyson 2008, in prep). 

The main LSST deliverables in the context of dark 
energy and matter will be 
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Fig. 24. — Simulations of the ratio of the measured galaxy power 
spectrum to a featureless reference power spectrum in various red- 
shift bins (shifted vertically for clarity), for the full LSST survey 
using very conservative photometric redshift errors. The several 
peaks visible in each curve are the signature of baryon acoustic 
oscillations (Eisenstein et al. 2005; Cole et al. 2005). LSST will 
measure the angular diameter distance over this redshift range with 
an accuracy of ~ 1% (Zhan & Knox 2006). 



• Multiple cosmic 2-point shear probes (Fig. 
There will be over 50 of these auto and cross power 
spectra constructed using 3-D shear tomography 
(correlation of shear of a galaxy at one redshift bin 
with another galaxy in a different redshift bin, av- 
eraged over all pairs of billions of galaxies; Jain & 
Taylor 2003). 

• Higher-order shear and galaxy statistics that can 
improve dark energy constraints and provide self- 
calibration of various systematics (Takada & Jain 
2004; Dolney, Jain & Takada 2006; Huterer et al 
2006). They are also probes of both primordial 
non-Gaussianities and those caused by non-linear 
structure. 

• Baryon Acoustic Oscillations in the galaxy angu- 
lar correlation functions. The standard ruler of the 
sound horizon at decoupling which is imprinted on 
the mass distribution at all redshifts provides a di- 
rect way of measuring the angular diameter dis- 
tance as a function of redshift (Fig. [Ml Eisenstein, 
Hu & Tegmark 1998; Cooray et al. 2001; Blake 
& Glazebrook 2003; Hu & Haiman 2003; Linder 
2003; Seo & Eisenstein 2003). LSST photo-z BAO 
will achieve percent level precision on the angular 
diameter distance at ~10 redshifts logarithmically 
spaced between z = 0.4 to 3.6 with this CMB- 
calibrated standard ruler (Fig. [^5)1 . When com- 
bined with CMB and weak lensing (WL) cosmic 
shear, this probe yields tight constraints on the 
dynamical behavior of dark energy. In particular, 
high-redshift BAO data can break the degeneracy 



between curvature and dark energy, and constrain 
ri/c to within 0.001, which is ten times better than 
the most accurate current result based on WMAP 
and SDSS data (Tegmark et al. 2006; Spergel et 
al. 2007). 

Supernovae. The two LSST samples are comple- 
mentary: the main survey will obtain light curves 
in 6 bands and photometric redshifts of about mil- 
lion Type la supernovae per year (permitting a 
search for a "third parameter" which, if not un- 
derstood, might introduce systematic error in lu- 
minosity evolution; Wood-Vasey et al. 2007), and 
the rapid sampling "mini-survey" of selected areas 
will yield well sampled light curves of tens of thou- 
sands of supernovae to a limiting redshift beyond 
one (leading to an independent test of dark energy 
dynamics; Riess et al. 2007). 

The distribution of strongest WL shear peaks with 
redshift. This is a simultaneous probe of the uni- 
versal mass function and the growth of structure, 
and provides an important constraint on dark en- 
ergy (Wang et al. 2005). LSST will find over 
200,000 such shear peaks on galaxy cluster scales. 

Unique constraints on anisotropy of cosmological 
parameters over the sky using SN, lensing clus- 
ters and lensing galaxies, which will test funda- 
mental cosmological assumptions of homogeneity 
and isotropy. 

The shape of the power spectrum of dark matter 
fluctuations measured by the LSST weak lensing 
maps will constrain the sum of neutrino masses 
with an accuracy of 0.04 eV or better (Cooray 1999; 
Song & Knox 2004; Hannestad, Tu & Wong 2006). 
The current best limit, as derived from a combi- 
nation of WMAP, the SDSS power spectrum, SN, 
and the Lya forest, is 0.17 eV (Seljak et al. 2006). 

Several millions of galaxy-galaxy lenses will pro- 
vide the needed statistics to constrain the dark 
matter mass function on small scales, and to probe 
dark matter halo profiles and substructure (Man- 
delbaum et al. 2006). 

LSST will discover several hundred strong lensing 
clusters in which background objects are multiply 
lensed. Cluster mass reconstruction based on the 
multiple image data and time delays between them 
(from SN or AGN activity) can probe the cluster in- 
ner mass profile, and can provide a separate test of 
cosmology, especially in cases with strongly lensed 
background objects at different redshift (Porciani 
& Madau 2000; Oguri & Kawano 2003). 

Time delays for QSOs multiply lensed by clusters 
as a function of redshift are an independent test 
of dark energy (Kundic et al. 1997). The natural 
timescale (many months to years) is well matched 
to the LSST survey. 



4.2. Taking an Inventory of the Solar System 

The small bodies of the Solar System, such as main- 
belt asteroids, the Trojan populations of the giant plan- 
ets and the Kuiper Belt objects, offer a unique insight 
into its early stages because they provide samples of the 
original solid materials of the solar nebula. Understand- 
ing these populations, both physically and in their num- 
ber and size distribution, is a key element in testing var- 
ious theories of the Solar System formation and the evo- 
lution. 

The baseline LSST cadence will result in orbital pa- 
rameters for several million moving objects; these will 
be dominated by main-belt asteroids, with light curves 
and multi-color photometry for a substantial fraction of 
detected objects. This dataset will yield 10 to 100 times 
more objects than are currently available with orbits, col- 
ors, and variability information. As already discussed in 
detail (§ l3.3.2p . LSST is capable of reaching the Congres- 
sional target completeness of 90% for PHAs larger than 
140 m, and will detect over 30,000 TNOs brighter than 
r ~ 24.5 using its baseline cadence. LSST will be capable 
of detecting objects like Sedna to beyond 100 AU, thus 
enabling in situ exploration far beyond the edge of the 
Kuiper belt at ~50 AU. Because most of these objects 
will be observed several hundred times, accurate orbital 
elements, colors, and variability information will also be 
available. 

Among the Solar System studies, these data will en- 
able: 

• Studies of the distribution of orbital elements for 
several million main-belt asteroids as a function 
of color-based taxonomy (see Fig. [26|) and size; 
size distributions of asteroid families (Parker et al. 
2008) and their correlations with age (Jedicke et 
al. 2004; Nesvorny et al. 2005); dynamical effects 
(Bottke et al. 2001); and studies of object shapes 
and structure using light curves (Pravec & Harris 
2000). 

• Studies of the distribution of orbital elements for 
about 100,000 NEOs as a function of color and 
size (Rabinowitz 1993; Dandy, Fitzsimmons & 
Collander-Brown 2002). Correlations with the 
analogous distributions for main-belt objects, and 
studies of object shapes and structure using light 
curves. 

• Studies of the distribution of orbital elements for 
about 100,000 Jovian Trojan asteroids as a func- 
tion of color and size (Jewitt, Trujillo & Luu 2000; 
Yoshida & Nakamura 2005; Szabo et al. 2007); the 
search for dynamical families (Knezevic & Milani 
2005); studies of shapes and structure using light 
curves. 

• Studies of the distribution of orbital elements for 
over 30,000 TNOs (see Fig. [27l) as a function of 
color and size; the search for dynamical families; 
studies of shapes and structure using light curves 
(Duncan & Levison 1995; Trujillo, Jewitt & Luu 
2001; Gladman et al. 2001; Bernstein et al. 2004; 
Elliot et al. 2005; Jones et al. 2006; Doressoundi- 
ram et al. 2007). 
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Fig. 26. — An example of color-based asteroid taxonomy. The 
figure shows the distribution of asteroids in the proper semi-major 
axis vs. sin{i) plane for 45,000 asteroids with colors measured by 
SDSS (Parker et al. 2008). The color of each dot is representative 
of the object's color. Note the strong correlation between aster- 
oid families (objects in distinct regions of the plane) and colors. 
There are at least five different taxonomic types distinguishable 
with SDSS measurements, and LSST color measurements of as- 
teroids will be more than twice as accurate (and will increase the 
number of objects by ~2 orders of magnitude). 

• An unbiased and complete census of both Jupiter- 
family and Oort-cloud comets; 6-band sub- 
arcsecond spatial profiles to a faint surface bright- 
ness limit; temporally resolved activity (Lowry et 
al. 1999; A'Hearn 2004). 

• The search for objects with perihelia at several 
hundred AU. For example, an object like Sedna 
(Brown, Trujillo & Rabinowitz 2004) will be de- 
tectable at 130 AU. 

• Studies of the solar wind propagation speed us- 
ing induced activity in a large sample of comets 
at different heliocentric distances (B. Vrsnak, priv. 
comm.) 



4.3. Exploring the Transient Optical Sky 

Time domain science will greatly benefit from LSST's 
unique capability to simultaneously provide large area 
coverage, dense temporal coverage, accurate color in- 
formation, good image quality, and rapid data reduc- 
tion and classification. Since LSST extends time-volume 
space a thousand times over current surveys (e.g., Mo- 
rokuma et al. 2008), the most interesting science may 
well be the discovery of new classes of objects. There are 
many projects that LSST data products will enable: 

• Detection and measurement of gamma-ray burst 
afterglows and transients (e.g., Zhang & Meszaros 
2004) to high redshift (-7.5). 

• Studies of optical bursters (those varying faster 
than 1 mag/hr) to r — 25 mag. 

• Gravitational micro-lensing in the Local Group and 
beyond (de Jong, Kuijken & Heraudeau 2008). 



10' 



S 10' 



E 

TO 



— I — I — I I 1 1 1 1 1 

5% albedo 
20% albedo 
50% aibedo 




Detection Distance at Opposition (AU) 

Fig. 27. — The LSST detection limits for distant solar system 
objects. Moving objects with diameters as small as 100m in the 
main asteroid belt and 100km in the Kuiper Be lt (TN Os) will be 
detected. Specialized deeper observations (see § 13. 2. 21 1 will detect 
TNOs as small as 10 km. Adapted from Jones ot al. (2007). 
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Fig. 28. — The phase space of cosmic explosive and eruptive 
transients as represented by their absolute V band peak brightness 
(j/ axis) and the event timescale (adapted from Kulkarni et al. 
2006). LSST will open up large regions of this phase space for 
systematic exploration by extending time- volume space 1000 times 
over existing surveys. 

• Studies of SN populations and parametrization of 
light curves (e.g., Hoeflich, Wheeler & Thielemann 
1998; Howeh et al. 2007). 

• Studies of dwarf novae, including their use as 
probes of stellar populations and structure in the 
Local Group (NeiU & Shara 2005). 

• A deep search for new populations of novae and 
supernova progenitors (Kulkarni et al. 2006; see 
FigEHl). 

• Search for stellar tidal disruptions by nuclear su- 
permassive black holes (Evans & Kochanek 1989; 
Gezari et al. 2008). 
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• A study of quasar variability using 1-2% accu- 
rate, multicolor light curves for ^2 million quasars, 
leading to constraints on the accretion physics (de 
Vries, Becker & White 2003; Vanden Berk et al. 
2004). 

• Optical identifications for transients detected at 
other wavelengths, from gamma-rays to radio, as 
well as by means other than electro-magnetic radi- 
ation (e.g., LIGO^o, LISA^i, ICECUBE^^) 

• The superb continuum light curves of AGN 
will enable economical "piggyback" reverberation- 
mapping efforts using emission lines (Liu et al. 
2008). These results will greatly broaden the 
luminosity-redshift plane of reverberation-mapped 
AGNs with black-hole mass estimates. For LSST 
data alone, the inter-band continuum lags will pro- 
vide useful structural information. 

• LSST data will allow good constraints on AGN life- 
times due to large sample size and survey lifetime 
(e.g.. Martini & Schneider 2003). 

4.4. Mapping the Milky Way 

The LSST will map the Galaxy in unprecedented de- 
tail, and by doing so revolutionize the fields of Galactic 
Astronomy and Near-field Cosmology. The great detail 
with which the Milky Way can be studied complements 
the statistical power of extra-galactic observations. The 
overarching goal of near-field cosmology is to use spatial, 
kinematic, and chemical data sets of stars to reveal the 
structure and evolution history of the Milky Way and its 
environment. LSST will reveal this fossil record in great 
detail and provide a Rosetta Stone for extragalactic as- 
tronomy by setting the context within which we interpret 
these broader data sets. Moreover, by providing a cen- 
sus of faint satellites and stellar streams in the halo, this 
map will offer a unique means to constrain the particle 
nature of dark matter because candidate supersymmet- 
ric particle dark matter models predict different mass 
clustering on small scales, with a corresponding range of 
mass profiles in low-mass systems. 

The LSST will produce a massive and exquisitely ac- 
curate photometric and astrometric data set for about 
10 billion Milky Way stars. The coverage of the Galactic 
plane will yield data for numerous star-forming regions, 
and the y band data will penetrate through the inter- 
stellar dust layer. Photometric metallicity measurements 
(see Figs. and [50)1 will be available for about 200 mil- 
lion main-sequence F/G stars which will sample the halo 
to distances of 100 kpc (Ivezic et al. 2008). No other ex- 
isting or planned survey will provide such a massive and 
powerful dataset to study the outer halo, including Gaia 
which is flux limited at r = 20, the SkyMapper which 
is flux limited to r = 22.6 (but will nicely complement 
SDSS in the Southern hemisphere), and Pan-STARRS 
which will not have the u band (which is necessary for 
estimating metallicity). The LSST in its standard sur- 
veying mode will be able to detect RR Lyrae variables 

http://www.ligo.caltech.edu 
http: / /lisa. nasa.gov 
http://icecube.wisc.edu 



(pulsating stars and standard candles) and classical no- 
vae (exploding stars and standard candles) at a distance 
of 400 kpc and hence explore the extent and structure of 
our own halo out to half the distance to the Andromeda 
galaxy. Thus, the LSST will enable studies of the distri- 
bution of main-sequence stars beyond the presumed edge 
of the Galaxy's halo (see Fig.[5T|). of their metallicity dis- 
tribution throughout most of the halo, and of their kine- 
matics beyond the thick disk/halo boundary. It will also 
obtain direct distance measurements via trigonometric 
parallax below the hydrogen-burning limit for a repre- 
sentative thin-disk sample. 

In addition to the study of hydrogen burning stars, 
LSST will uncover the largest sample of stellar remnants 
to date. Over 97% of all stars eventually exhaust their 
fuel and cool to become white dwarfs. Given the age 
of the Galactic halo, most of the mass in this compo- 
nent is now tied up in these remnant stars (e.g., Alcock 
et al. 2000) and therefore their discovery directly con- 
straints the Galactic mass budget. These large popula- 
tions of disk and halo white dwarfs will provide unprece- 
dented constraints on the luminosity function of these 
stars, which will directly yield independent ages for the 
Galactic disk and halo (e.g., through the initial-final mass 
relation, Kalirai et al. 2008). 

The spatial coverage of LSST naturally targets both 
field stars and star clusters. To date, no systematic sur- 
vey of the stellar populations of Southern hemisphere 
clusters has been carried out (e.g., such as the CFHT 
Open Star Cluster Survey or the WIYN Open Star Clus- 
ter Survey in the North; Kalirai et al. 2001; Mathieu 
2000). Multiband imaging of these co-eval, co-spatial, 
and iso-metallic systems will provide vital insights into 
fundamental stellar evolution. For example, the depth 
of LSST will enable us to construct luminosity and mass 
functions for nearby open clusters down to the hydrogen 
burning limit. Variations in the initial mass function will 
be studied as a function of environment (e.g., age and 
metallicity). The wide-field coverage will also allow us 
to track how the stellar populations in each cluster vary 
as a function of radius, from the core to beyond the tidal 
radius. Fainter remnant white dwarfs will be uncovered 
in both open and globular clusters (the nearest of which 
are all in the South) thereby providing a crucial link to 
the properties of the now evolved stars in each system. 

In summary, the LSST data will revolutionize studies 
of the Milky Way and the entire Local Group (Saha et al. 
2007). We list a few specific Galactic science programs 
that LSST wiU enable: 

• High-resolution studies of the distribution of stars 
in the outer halo in the six-dimensional space 
spanned by position, metallicity and proper mo- 
tions (e.g., Girard et al. 2006; BeU et al. 2007; 
Juric et al. 2008). 

• The faintest ever search for halo streams, and 
galaxy satellites and intergalactic stars over much 
of the Local Group (e.g., Grillmair 2006ab; Ibata 
et al. 2007; Walsh, Jerjen & Willman 2007; Be- 
lokurov et al. 2007a). 

• Deep and highly accurate color-magnitude dia- 
grams for over half of the known globular clusters, 
including tangential velocities from proper motion 
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Fig. 29. — The g — r vs. u — g color-color diagram for two mil- 
lion stars measured by SDSS (adapted from Smolcic et al. 2004). 
Accurate multi-color photometry contains information that can be 
used for source classification and determination of detailed stellar 
properties such as effective temperature and metallicity. LSST will 
enable such measurements for billions of stars. 
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Fig. 30. — The median metallicity map for 2.5 million main- 
sequence F-type stars within 10 kpc from the Sun (adapted from 
Ivezic et al. 2008). The metallicity is estimated using u — g and 
g—r colors measured by SDSS. The position and size of the mapped 
region, relative to the rest of the Milky Way, is illustrated in the top 
right corner, where the same map is scaled and overlaid on an image 
of the Andromeda galaxy. The gradient of the median metallicity 
is essentially parallel to the Z axis, except in the Monoceros stream 
region, as marked. LSST will extend this map out to 100 kpc, using 
a sample of over 100 million main-sequence F stars. 
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Fig. 31. — A simulation of the outer regions of the Milky Way 
(left), compared to the current state-of-the-art data (right), shown 
on the same spatial scale. The data panel shows the number density 
multiplied by the cube of the Galactocentric radius (logarithmic 
scale with dynamic range of 1000, from blue to red), for ~1000 
SDSS RR Lyrae stars within 10 degrees of the Sgr dwarf tidal 
stream plane (Ivezic et al. 2003). The same color coding was used 
to visualize stellar number density for a Milky Way type galaxy 
simulation from Bullock and Johnston (2005), shown on the left. 
Set with an ACDM merger history, these simulations track the 
accretion and disruption of hundreds of dwarf galaxies into Milky- 
Way size halos. With LSST, RR Lyrae stars will be found beyond 
the presumed Milky Way tidal radius (~300 kpc, white circle), and 
much more numerous main-sequence stars will trace the structure 
out to 100 kpc (dashed circle), with similar fidelity as shown in 
Figure lSOl The latter distance range can now be probed only using 
RR Lyrae stars and other rare non-main-sequence stars. 

measurements (Casetti-Dinescu et al. 2007; Clem, 
VandenBerg & Stetson 2008). 

• Mapping the metallicity, kinematics and spatial 
profile of the Sgr dwarf tidal stream (e.g., Ibata 
et al. 2001; Majewski et al. 2003; Law, Johnston 
& Majewski 2005). 

• Studies of the dumpiness of the gravitational po- 
tential in the Galaxy using fragile wide-angle bi- 
naries selected with the aid of trigonometric and 
photometric parallaxes, and common proper mo- 
tion (e.g., Yoo, Chaname & Gould 2004). 

• Detailed studies of variable star populations; 2% or 
better accurate multicolor light curves will be avail- 
able for a sample of at least 50 million variable stars 
(Sesar et al. 2007), enabling studies of cataclysmic 
variables, eclipsing binary systems, rare types of 
variables, etc. 

• Discovery of rare and faint high proper motion ob- 
jects: probing the faint end of the stellar mass func- 
tion (Lepine 2008), and searching for free-floating 
planet candidates (Lucas & Roche 2000). 

• Direct measurement of the faint end of the stellar 
luminosity function using trigonometric parallaxes 
(Gratton et al. 1997). A complete census of the 
solar neighborhood to a distance of 100 pc based on 
trigonometric parallax measurements for objects as 
faint as = 17. For example, LSST will deliver 
10% or better distances for a sample of about 2,500 
stars with 18< Mr <19. There are only a handful 
of such stars known today (and Gaia will detect 
fewer than 100). 
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Galaxy Environment 




Fig. 32. — A comparison of the distribution of galaxies in 
luminosity— color— density space measured by SDSS (left) and a 
model based on the Millennium simulation (right). The linearly- 
spaced contours outline the distribution of volume-limited sample 
of galaxies in the plotted diagram, and the color-coded background 
shows the median environmental density (computed using the ten 
nearest neighbors) for galaxies with the corresponding luminosity 
and color. Such multi-variate distributions encode rich information 
about formation and evolution of galaxies. Galaxies detected by 
SDSS are representative of the low-redshift Universe (the median 
redshift is ~0.1). The LSST will enable such studies as a function 
of redshift, to z ~2. Adapted from Cowan & Ivczic (2008). 

• Using their rich molecular band structure, the 
LSST colors (especially the z — y color) and faint 
limit will permit the separation of halo M sub- 
dwarfs from disk M dwarfs (West et al. 2004; Lep- 
ine & Scholz 2008). 

• Planetary transits: the data set may include a mini 
survey of 600 sq.deg. of sky which will collect about 
40 hour-long sequences of 200 observations each 
over a 4-month period (see § I3.2.2P . There would 
be an additional 800 observations over 10 years of 
these same fields as part of the main survey. With 
about 10 million or more stars in the sample (de- 
pending on where the fields are placed) , this would 
provide an excellent material to study the planet 
frequency as a function of stellar type, metallicity, 
and distance from the Galactic plane (e.g., Char- 
bonneau et al. 2000; Konacki et al. 2003). 

• A complete census of AGB stars in the Galaxy by 
searching for resolved envelopes and optical identi- 
fications of IR counterparts (e.g., from IRAS and 
Spitzer surveys), and by using long-term variability 
and color selection (Ivezic 2007c). 

• A complete census of faint populations in nearby 
star forming regions using color and variability se- 
lection (e.g., Briceno et al. 2005). 

• High-resolution three-dimensional studies of inter- 
stellar dust using 5-color SEDs of main sequence 
stars (McGehee 2004; Meyer et al. 2005). 



4.5. Additional Science Projects 

The experience with any large survey (e.g., SDSS, 
2MASS, GALEX, to name but a few) is that the most 
interesting science that comes out is often unrelated to 
the main science drivers, and is unanticipated at the time 
the survey is designed. LSST will enable even more di- 
verse science than encompassed by the four themes that 
drive the system design. We list a few anticipated major 
programs. 
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Fig. 33.— The LSST will deliver AGN sky densities of 1000- 
4000 deg~^, which is competitive with the deepest pencil-beam 
surveys to date, such as the Chandra Deep Fields (top panel). The 
total LSST AGN yield, selected using colors and variability, should 
be 20-80 million objects. The bottom panel shows the expected 
distribution of these objects in the absolute magnitude vs. redshift 
plane, color-coded by the probability for an object to be detected 
as variable after 1 year of observations. Note that quasars will 
be detected to their formal luminosity cutoff (A/ < —23) even at 
redshifts of ~5. Adapted from Brandt et al. (2007). 

• Detailed studies of galaxy formation and evolu- 
tion using their distribution in luminosity-color- 
morphology space as a function of redshift. For ex- 
ample, LSST will enable studies of the rest-frame 
UV emission, similar to those based on GALEX 
data for local galaxies, to a redshift of ~2. These 
studies project onto many axes (see Fig. [5^ : 

— the evolution of the galaxy luminosity func- 
tion with redshift, as a function of morphology 
and color 

— the evolution of the galaxy color distribution 
over a wide range of rest-frame wavelengths, 
and as a function of luminosity and morphol- 
ogy 

— bulge-disk decomposition, as a function of lu- 
minosity and color, over a large redshift range 
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— detailed distribution of satellite galaxies in 
luminosity-color-morphology space as a func- 
tion of luminosity, color, and morphology of 
the primary galaxy 

— correlations of luminosity, color and morphol- 
ogy with local environment, and as a function 
of redshift 

— the properties of galaxy groups as a function 
of cosmic time. 

• AGN census to a very faint luminosity limit and 
to a large redshift limit, yielding 20-80 million ob- 
jects (see Fig. [33l). By pushing substantially fur- 
ther down the AGN luminosity function over a very 
large solid angle, LSST data will lead to a much 
clearer understanding of black-hole growth during 
the first Gyr. For example, LSST should discover 
~1000 AGNs at z - 6.5 - 7.5 at the end of the 
dark ages, representing a dramatic increase over 
the present samples (Brandt et al. 2007). 

• The first half-sky survey of ultra low surface bright- 
ness objects, with photometric redshift informa- 
tion. The currently available samples are highly 
incomplete, especially in the Southern Hemisphere 
(see Fig. 7 in Belokurov et al. 2007a). 

• Search for strong gravitational lenses to a faint sur- 
face brightness limit (e.g., Bartelmann et al. 1998; 
Tyson, Kochanski & DeU'Antonio 1998; Belokurov 
et al. 2007b). 

• Multi-wavelength studies of faint optical sources 
using X-ray, IR and radio data. For example, 
SDSS detected only 1/3 of all 20cm FIRST sources 
(Becker et al. 1995) because it was too shallow by 
~4 mag for a complete optical identification. 

4.5.1. Synergy with other projects 

LSST will not operate in isolation and will greatly ben- 
efit from other precursor and coeval data at a variety of 
wavelengths, depths, and timescales. For example, most 
of the Celestial Sphere will be covered to a limit sev- 
eral magnitudes fainter than LSST saturation (r ~ 16), 
first by the combination of SDSS, PSl and SkyMapper 
(Keller et al. 2001), and then by the Gaia survey. The 
SkyMapper survey will obtain imaging data in the south- 
ern sky to similar depths as SDSS, the Pan-STARRS sur- 
veys will provide multi-epoch data deeper than SDSS in 
the northern sky, and the Dark Energy Survey (Flaugher 
2006) will scan ^5000 deg^ to similar depth as Pan- 
STARRS in the southern sky. Despite the lack of the 
u band data and a shallower depth by about a magni- 
tude, the Pan-STARRS surveys will represent a valuable 
complement to LSST in providing an all-sky coverage to 
a limit fainter than that of SDSS and SkyMapper. LSST 
and Gaia will be complementary datasets for studying 
the Milky Way in the multi-dimensional space of three- 
dimensional positions, proper motions and metallicity. 
The Gaia survey will provide calibration checks at the 
bright end for proper motions and trigonometric paral- 
lax measurements by LSST, and LSST will extend the 
Gaia survey four magnitudes deeper. The LSST data 



stream will invigorate subsequent investigations by nu- 
merous other telescopes that will provide additional tem- 
poral, spectral and spatial resolution coverage. 

LSST will provide a crucial complementary capabil- 
ity to planned space experiments operating in other 
bands, such as the Gamma-Ray Large Area Space Tele- 
scope (GLAST^^) and the Black Hole Finder Einstein 
Probe Mission. The Laser Interferometer Space Antenna 
(LISA) is designed to study black hole mergers through 
the gravitational wave outburst that should be emitted 
during the final stages of such events. LSST has the 
potential to measure the electromagnetic signal that ac- 
companies the gravitational wave emission, thereby pro- 
viding a definite position on the sky for the system, which 
is crucial for follow-up observations. LSST will also make 
the archives for billion-dollar class space missions such as 
Chandra, XMM-Newton, Spitzer, etc., much more valu- 
able, because deep optical multi-color data will enable 
massive photometric follow-up studies of sources from 
these missions. Last but not least, the huge samples 
of various astronomical source populations will yield ex- 
tremely rare objects for follow-up by powerful facilities 
such as JWST (Gardner et al. 2006) and GSMT (Kan 
& Antebi 2003). 

In summary, the diversity of science enabled by LSST 
will be unparalleled, extending from the physics of grav- 
ity and the early Universe to the properties of "killer" 
asteroids (PHAs). While there are other projects that 
aim to address some of the same science goals, no other 
project matches this diversity and LSST's potential im- 
pact on society in general. 

5. COMMUNITY INVOLVEMENT 

LSST has been conceived as a public facility: the 
database that it will produce, and the associated object 
catalogs that are generated from that database, will be 
made available to the scientific community and to the 
public at large with no proprietary period. The LSST 
data management system will provide user-friendly tools 
to access this database and to support user-initiated 
queries, run on LSST computers, either at the archive 
facility or at the data access centers. We expect that 
many, perhaps even the majority of LSST discoveries will 
come from research astronomers with no formal affilia- 
tion to the project, from students, and from interested 
amateurs, intrigued by the accessibility to the Universe 
that this facility uniquely provides. For example, because 
of its interest in the large public data access of the LSST 
program, the Google Corporation has recently joined the 
LSST team. 

The SDSS provides a good example for how the scien- 
tific community can be effective in working with large, 
pubhcly available astronomical data sets. The SDSS 
has put out a series of large incremental data releases 
via a sophisticated database, roughly once per year, to- 
gether with a paper describing the content of each data 
release, and extensive on-line documentation giving in- 
stru ctions on down loading the catalogs and image data 
(see |http://www. sdss.org). Over 60% of the more than 
1800 refereed papers based on SDSS data to date (May 
2008) have been written by people from outside the 
project, and include many of the most high-profile re- 

^■^ http://glast.gsfc.nasa.gov 
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suits that have come from the facihty. 

Nevertheless, it is clear that many of the highest pri- 
ority LSST science investigations will require organized 
teams of professionals working together to optimize sci- 
ence analyses and to assess the importance of systematic 
uncertainties on the derived results. To meet this need, 
ten science collaborations have already been established 
by the project in core science areas (Strauss 2007). 

As of the time of this contribution, there are close 
to 200 participants in these collaborations, mostly from 
LSST member institutions. As the project proceeds into 
development, an open call for applications for member- 
ship will be issued to the professional community at large. 
It is anticipated that the LSST science collaborations will 
perform their analyses using LSST computational facili- 
ties, and the system has been sized accordingly. 

As we design our observing strategies, we are actively 
seeking and implementing input by the wide community 
of scientists in the LSST Science Collaborations. Their 
task, among others, is to enhance the LSST science case 
and provide advice on how to optimize their science with 
choices in cadence, software, and data systems. The Sci- 
ence Collaborations are expected to play a major role in 
the system optimization during the commissioning pe- 
riod. We will continue this review when the survey starts, 
and will maintain our design of flexible cadence structure. 

6. EDUCATIONAL AND SOCIETAL IMPACTS 

The impact - and enduring societal significance - of the 
LSST will exceed its direct contributions to advances in 
physics and astronomy. With an open database available 
in real time to the non-specialist public, LSST has the 
capacity to literally bring the Universe home to everyone. 
The "movie of the Universe" will offer a new window on 
the sky for curious minds of all ages. Any interested 
amateur, anywhere in the world, will be able to log on 
to his or her computer and look at the latest images of 
any part of the sky^''. In the context of open source de- 
velopment, as we have witnessed recently in community 
contributions to Google Sky, the value added for public 
interest and interaction is enormous. LSST will allow 
people to see that they are part of a dynamic Universe. 
User-friendly visualization interfaces will enable new re- 
search as well as being a vehicle for education and public 
involvement. Elementary school students in Earth sci- 
ence classes can participate in the discovery of new as- 
teroids. High school and college students can discover su- 
pernovae and measure their light curves, thereby learning 
about both measurement processes and cosmology. The 
implications for all levels of scientific, engineering, and 
technical education cannot be overstated. 

The LSST Education and Public Outreach (EPO) pro- 
gram will provide tools for the public, museum visitors, 
and students not only to examine and understand the 
achievements of LSST, but also to participate in the pro- 
cess of scientific discovery. With partners and collab- 
orators such as the Google Corporation, the American 
Museum of Natural History, the Las Cumbres Observa- 
tory Global Telescope Network, and others, a number 

At the angular resolution of LSST, it would take over 3 million 
HDTV sets (with 1080x1920 pixel resolution) to display a single- 
epoch LSST image of half the Celestial Sphere (and about 1000 
such images will be obtained over 10 years). 



of innovative programs have been devised: an internet- 
based public outreach system called involving citizen- 
scientists, advanced visualization demonstrations at the 
Hayden and Adler planetaria, formal educational pro- 
grams focused on solar system science and cosmology, 
teacher training, and the development of new software 
EPO tools. LSST is perfectly positioned to bring such a 
transformative EPO program to fruition. 

7. SUMMARY AND CONCLUSIONS 

Until recently, most astronomical investigations have 
focused on small samples of cosmic sources or individ- 
ual objects. Over the past decade, however, advances in 
technology have made it possible to move beyond the tra- 
ditional observational paradigm and to undertake large- 
scale sky surveys, such as SDSS, 2MASS, GALEX and 
many others. This observational progress, based on a 
synergy of advances in telescope construction, detectors, 
and above all, information technology, has had a dra- 
matic impact on nearly all fields of astronomy, many ar- 
eas of fundamental physics, and society in general. The 
LSST builds on the experience of these surveys and ad- 
dresses the broad goals stated in several nationally en- 
dorsed reports by the U.S. National Academy of Sciences. 

The realization of the LSST involves extraordinary en- 
gineering and technological challenges: the fabrication 
of large, high-precision optics; construction of a huge, 
highly-integrated array of sensitive, wide-band imaging 
sensors; and the operation of a data management facility 
handling tens of terabytes of data each day. The design 
and development effort, which has been underway since 
2006 and will continue through the onset of construction, 
includes structural, thermal, and optical analyses of all 
key hardware subsystems, vendor interactions to deter- 
mine manufacturability, explicit prototyping of high-risk 
elements, prototyping and development of data manage- 
ment systems, and extensive systems engineering studies. 
Over 100 technical personnel at a range of institutions are 
currently engaged in this program. 

In September 2007, LSST passed the NSF Conceptual 
Design Review for construction, which puts it on track 
for the Preliminary Design Review around the end of 
2008. The casting of the primary/tertiary mirror is in 
progress^^ and it reached a peak temperature of approx- 
imately 1165 °C on March 29'^, 2008. The LSST design 
and development phase will continue for another 3 years 
under NSF support along with in-kind contributions and 
private gifts. It is proposed that the U.S. Department 
of Energy support the cost of constructing the camera. 
LSST has high visibility in the high-energy physics com- 
munity, both at universities and government laborato- 
ries. The project is scheduled to have first light in 2014 
and the beginning of survey operations in 2015. 

The LSST survey will open a movie-like window on 
objects that change brightness, or move, on timescales 
ranging from 10 seconds to 10 years. The survey will 
have a data rate of about 30 TB/night (more than one 
complete Sloan Digital Sky Survey per night), and will 
collect over 60 PB of raw data over its lifetime, resulting 
in an incredibly rich and extensive public archive that 
will be a treasure trove for breakthroughs in many ar- 
eas of astronomy and physics. About 10 billion galaxies 

Please see |http: //www.lsst.org/News/highfire-event.shtml| 
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and a similar number of stars will be detected - for the 
first time in history, the number of cataloged celestial ob- 
jects will exceed the number of living people! About a 
thousand observations of each position across half of the 
Celestial Sphere will represent the greatest movie of all 
time: assuming HDTV resolution (2.1 Mpix/frame) with 
30 frames per second, and 270 color images per position 
(e.g., the ugr and izy color-composites, c.f. Table 1), it 
would take about eleven uninterrupted months to view 
this color movie (equivalent to about 4,000 feature films). 

The cost estimate for LSST is $389M for construction 
and $37M per year for operations (in 2006 U.S. dollars). 
LSST will be much more efficient than the current state- 
of-the-art massive optical survey, the SDSS. For exam- 
ple, for each $100 spent, the SDSS obtained single-epoch 
ugriz photometry for '^200 sources and a spectrum for 
a single source. In contrast, for $100 LSST will obtain 
ugrizy photometry for more than 2,000 sources (to a 
limit 5 magnitudes deeper than SDSS) and a 1000-point 
light curve for ^--^200 sources. Compared to the SDSS, 
LSST will produce about 2,000 times larger data volume 
with an increase in total cost of less than a factor of ten. 

LSST has been conceived as a pubhc facility: the 
database that it will produce, and the associated object 
catalogs that are generated from that database, will be 
made available to the world's scientific community and 
to the public at large with no proprietary period. The 
software which creates the LSST database will be open 
source. LSST will be a significant milestone in the glob- 
alization of the information revolution. LSST will put 
terabytes of data each night into the hands of anyone 
that wants to explore it, and in some sense will become 
an internet telescope: the ultimate network peripheral 
device to explore the Universe, and a shared resource for 
all humanity. 
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